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NUMERICAL  AND  ALGEBRAIC 


Gain  in  decibels  of  P2  relative  to  Pi 


G  =  101og10(P2/Pi). 


To  within  two  percent 

(2? r)1/2  «  2.5;  7 r2  «  10;  e3  »  20;  210  »  103. 


Euler-Mascheroni  constant1  7  =  0.57722 


Gamma  Function  T(a;  +  1)  =  ®r(a;): 


r(i/6) 

r(l/5) 

r(i/4) 

r(i/3) 

F(2/5) 

r(i/2) 


5.5663 

4.5908 

3.6256 

2.6789 

2.2182 

1.7725  =  y/TV 


r(3/5)  =  1.4892 
r(2/3)  =  1.3541 
r(3/4)  =  1.2254 
r(4/5)  =  1.1642 
r(5/6)  =  1.1288 
r(l)  =  1.0 


Binomial  Theorem  (good  for  |  x  |<  1  or  a  =  positive  integer): 


(1  +  X)«  =  ^2  Q*k  =  1  +  ax  +  Q(a2!  ■— x2  + 


a(a  —  1)  2  .  a  (a  —  l)(a  -  2)  3 


3! 


k—0 


Rothe-Hagen  identity2  (good  for  all  complex  x,  y,  z  except  when  singular): 

71 

x  (X  +  kz\  y  /y  +  (n  —  k)z\ 


Ex  /x  +  kz\ 
x  +  kz  V  k  ' 


k= 0 


y  +  (n  —  k)z  V  n  —  k 

x  +  y  (X  +  y  + 


x  +  y  +  nz 


/x  +  y-t  nz\ 
V  n  ' 


Newberger’s  summation  formula3  [good  for  y,  nonintegral,  Re  (a  +  (3)  >  —  1]: 

oo 

(  1)  J OL  ~  771  (^)  J (3  +  ^71  ^ 


E 


n—  —  oo 


n  +  /i 


Sin  /Z7T 


Jql+ 7/x  —  7M  (^O* 
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VECTOR  IDENTITIES4 

Notation:  /,  g,  are  scalars;  A,  B,  etc.,  are  vectors;  T  is  a  tensor;  /  is  the  unit 
dyad. 

(1) ABxC  =  AxBC  =  BCxA  =  BxCA  =  CAxB  =  CxAB 

(2)  A  x  (B  x  C)  =  (C  x  B)  x  A  =  (A  •  C)B  -  (A  •  B)C 

(3)  A  x  (B  x  C)  +  B  x  (C  x  A)  +  C  x  (A  x  B)  =  0 

(4)  (A  x  B)  •  (C  x  D)  =  (A  •  C)(B  •  D)  -  (A  •  D)(B  •  C) 

(5)  (A  x  B)  x  (C  x  D)  =  (A  x  B  •  D)C  -  (A  x  B  •  C)D 

(6)  V(fg)  =  V(gf)  =  fVg  +  gVf 

(7)  V  •  (/A)  =  /V  •  A  +  A  •  V/ 

(8)  V  x  (/A)  =  /V  x  A  +  V/  x  A 

(9)  V-(AxB)  =  B  •  V  x  A-A'V  xB 

(10)  V  x  (A  x  B)  =  A(V  •  B)  -  B(V  •  A)  +  (B  •  V)A  -  (A  •  V)B 

(11)  A  x  (V  x  B)  =  (VB)  •  A  -  (A  •  V)B 

(12)  V(A  •  B)  =  A  x  (V  x  B)  +  B  x  (V  x  A)  +  (A  •  V)B  +  (B  •  V)A 

(13)  V2/  =  V- V/ 

(14)  V2A  =  V(V  •  A)  -  V  X  V  x  A 

(15)  V  x  V/  =  0 

(16)  V  •  V  x  A  =  0 

If  ei,  e2,  e3  are  orthonormal  unit  vectors,  a  second-order  tensor  T  can  be 
written  in  the  dyadic  form 

(17)  T  =  r«e(e, 

In  cartesian  coordinates  the  divergence  of  a  tensor  is  a  vector  with  components 

(18)  (V-T)i  =  '£.(dTji/dxj) 

[This  definition  is  required  for  consistency  with  Eq.  (29)].  In  general 

(19)  V  •  (AB)  =  (V  •  A)B  +  (A  •  V)B 

(20)  V  -  (/T)  =  V/-T+/V-T 
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Let  r  =  \x  +  jy  +  kz  be  the  radius  vector  of  magnitude  r,  from  the  origin  to 


the  point  x:  y,  z. 

Then 

(21) 

V  •  r  =  3 

(22) 

Vxr  =  0 

(23) 

Vr  =  r/r 

(24) 

V(l/r)  = 

-r/r3 

(25) 

V • (r/r3) 

=  4-7r<5(r) 

(26) 

Vr  =  / 

If  V  is  a  volume  enclosed  by  a  surface  S  and  dS  —  ndS,  where  n  is  the  unit 
normal  outward  from  V, 


(27) 

(28) 

(29) 

(30) 

(31) 

(32) 


dVV f  =  /  dSf 


L 


dVV  •  A  =  /  dS  ■  A 


I, 


l 
L 

/  dW-T  =  dS  T 
J  v  J  s 

L 
L 


dVV  x  A  =  /  dSxA 


dV(fV2g-gV2f)=  /  dS-(fVg-gVf) 


J, 


L 


dV(  A  -  VxVxB  —  B-VxVxA) 


l 


=  /  dS-(BxVxA-AxVxB) 


If  5  is  an  open  surface  bounded  by  the  contour  C,  of  which  the  line  element  is 
dl, 


(33) 


dS  x  V/  = 


dlf 
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(34) 

(35) 

(36) 


1, 

l 

L 


dS  ■  V  x  A  =  (b  dl  •  A 

c 

(dS  x  V)  x  A  =  j)  dl  x  A 
dS  ■  (V/  x  V<?)  =  (j)  fdg  =  - 


DIFFERENTIAL  OPERATORS  IN 
CURVILINEAR  COORDINATES5 


Cylindrical  Coordinates 

Divergence 


V  •  A  = 


Ar)  + 

r  or 


1  dA$ 

r  d(j) 


+ 


dAz 

dz 


Gradient 

(V/)r 


dr  ’ 


(V/)z 


dl 

dz 


Curl 


(V  x  A)r 


1  dAz 
r  d(j) 


dA $ 

~dT 


(V  x  A)* 


dAr  dAz 
dz  dr 


(V  x  A)* 


I±  (rxn-i  — 

r9rl  r  d(p 


Laplacian 

r  dr  \  dr  J  r2  d(f>2  dz2 
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Laplacian  of  a  vector 

(V2A)r  =  V2A.-4^!±-4 
rz  ocp  rz 


z  2  *  \  _ 2  ,,  2  9Ar  Aa) 

(V  A)*=v  At  +  -—~  — 


(V2A)*  =  v2a: 


Components  of  (A  •  V)B 

(A- VB  )r  =  Ar^  +  +  Az9Br  A ^ 


<9r  r  dp 


dz 


(A.VB)^Ar2pL  +  ±tB±  +  A^  +  A^ 


dr  r  dip 


dz 


Divergence  of  a  tensor 

—  T-N  _  1  ^  /  rp  \  |  1  Tfjxj, 

(V  •  /  )r  —  —  -jr-(rTrr)  +  + 

r  dr  r  dp  dz  r 


i  d  ,  _  .1  dT^  arz<^  tv 


(-V'T^  =  -rirPrTr*)  +  -r-W+  8 , 


+ 


i  a  .  m  .  i  axv  ar,* 

(V  •  T)z  —  -  —  (rTrz)  H - — — | — — — 

r  c/r  r  a<p  az 
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Spherical  Coordinates 


Divergence 


A  =  -4^-(r2Ar)  + 

r z  or 


d 


r  sin  6  dO 


—  (sin  6Ae)  + 


1  dA $ 

rsinO  d(f) 


Gradient 


(X7n  =  21.  mn  - -21.  mn  =  1  21 

(  f)r  Qr,  (  f)e  r  Qe ,  (  fU  rsin050 


Curl 


(V  x  A)r  =  — -2- (sin  OAtf,) - ?— r  -jrr 

r  sin  0  a0  r  sin  0  ocp 


(V  x  A)e  =  — —  ^  -  —  -2-(rA^) 

r  sin  0  o0  rar 


(V  x  A),/,  =  -^-(rAo)  -  ~-£jr 


Laplacian 


V2/ 


=  ±d_  (  *df\  1  d  (sin0d£\  _ 1 

r2  dr  \  dr  )  r2  sin  6  d6  \  dd  )  r 2  sii 


a2/ 

sin2  6  d<f> 2 


Laplacian  of  a  vector 

2  A  n.  2  4  2Ar  2  aA#  2  cot  2  aA^, 

(V  ==  \7  A>p 


r2  r2  a$ 


r2  sin  0  d<p 


(V2A)9  =  v2a»  +  A  _  _Af _ 2cos()  gd* 

r2  a#  r2  sin2  0  r2  sin2  6  d(f> 


(V2A)^>  =  v2a*  - 


A^ 


+ 


2  aAr  2  cos  0  aA# 


+ 


r2  sin2  0  r2  sin  6  d$>  r2  sin2  0  d(f> 
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Components  of  (A  •  V)B 


(A.  •  VB)r  —  Ar 


dBr 

dr 


Ag  dBr  A<j>  dBr 

+  r  d9  '  r  sin  0  d(j> 


AgBg  +  AtpBj, 

r 


(A-VB  )e  =  Ar^ 
or 

(A- VB )<p  =  Ar^ 
Divergence  of  a  tensor 


Aq  dB$  A $  8Bq  AqBt  cot OA^Bfj) 

r  86  rsin#  d<f)  +  r  r 

^  Ao  dB$  A $  dBcf)  A^Br  cot  QA^Bp 

^  r  dO  ^  r  sin  6  dcj)  r  r 


(V  *  T)r  —  —  —  (r  ^rr)  + 
rz  or 


d 


r  sin  0  dO 


(sin  OTor ) 


1  dTffrr  T$e  + 


rsin0  d(j) 


(V  •  T)e  = 


9  ( r2Tre )  + 


d 


dr 


r  sin  6  dO 


—  (sin  OTgo ) 


1  dT^g  Tgr  cot  OTfj, (f) 
r  sin  6  d<f>  r  r 


(  V7 


rsin#  d(f> 


dT '4,4,  T$r  +  cot  6T$e 
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DIMENSIONS  AND  UNITS 


To  get  the  value  of  a  quantity  in  Gaussian  units,  multiply  the  value  ex¬ 
pressed  in  SI  units  by  the  conversion  factor.  Multiples  of  3  in  the  conversion 
factors  result  from  approximating  the  speed  of  light  c  =  2.9979  X  1010  cm/sec 
«3x  1010  cm/sec. 


Dimensions 

Physical 

Sym- 

SI 

Conversion 

Gaussian 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

Capacitance 

C 

,2  2 
t  q 

Tfll^ 

l 

farad 

9  X  1011 

cm 

Charge 

Q 

q 

m}/2lz/2 

t 

coulomb 

3  X  109 

statcoulomb 

Charge 

density 

q 

mV2 

coulomb 

/m3 

3  x  103 

statcoulomb 

/cm3 

P 

i 3 

PPt 

Conductance 

tq2 
ml 2 

l 

t 

siemens 

9  x  1011 

cm/sec 

Conductivity 

a 

tq2 
ml 3 

1 

7 

siemens 

/m 

ampere 

9  x  109 

sec-1 

Current 

i,i 

<? 

m1/2/3/2 

3  x  109 

statampere 

t 

P 

Current 

density 

j,j 

q 

m1/2 

ampere 

/m2 

3  x  105 

statampere 

/cm2 

IH 

PPP 

Density 

p 

m 

p 

m 

~P 

kg/m3 

10-3 

g/cm3 

Displacement 

D 

q 

«^l/2 
m  ' 

coulomb 

/m2 

12t r  X  105 

statcoulomb 

/cm2 

P 

p/n 

Electric  field 

E 

ml 

m1/2 

volt/m 

-  x  10-4 

3 

statvolt/cm 

P~q 

P/H 

Electro- 
mot  ance 

e, 

Emf 

ml2 

m1/2l1/2 

volt 

-  x  10-2 

3 

statvolt 

t2q 

t 

Energy 

U,W 

ml2 

t2 

ml 2 

P 

joule 

107 

erg 

Energy 

density 

w,  e 

m 

Jp 

m 

Jt2 

joule/m3 

10 

erg/cm3 
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Dimensions 

Physical 

Sym- 

SI 

Conversion 

Gaussian 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

ml 

ml 

dyne 

Force 

F 

~w 

newton 

105 

Frequency 

/,  v 

i 

t 

i 

t 

hertz 

1 

hertz 

ml2 

t 

i  x  l<Tn 

9 

sec/cm 

Impedance 

Z 

tq 2 

1 

ohm 

Inductance 

L 

ml2 

q2 

t 2 
i 

henry 

1  _ ii 

-  x  10  11 

9 

sec2 /cm 

Length 

l 

l 

i 

meter  (m) 

102 

centimeter 

(cm) 

Magnetic 

H 

q 

m1/2 

ampere- 

4t r  X  10-3 

oersted 

it 

ll/2t 

intensity 

turn/m 

ml 2 

m1/2Z3/2 

108 

maxwell 

Magnetic  flux 

$ 

tq 

m 

t 

weber 

m1'2 

Magnetic 

induction 

n 

tesla 

104 

gauss 

JZ> 

tq 

l'/2t 

Magnetic 

m,  fi 

m1/2^2 

ampere-m2 

103 

oersted— 

4. 

moment 

t 

t 

3 

cm 

Magnetization 

M 

q 

m 1^2 

ampere- 

1— l 

O 

1 

CO 

oersted 

it 

turn/m 

Magneto- 

M, 

q 

m1'2^'2 

ampere- 

47T 

gilbert 

t2 

10 

motance 

Mmf 

t 

turn 

Mass 

m,  M 

m 

m 

kilogram 

103 

gram  (g) 

ml 

t 

ml 

T 

(kg) 

Momentum 

P,P 

kg-m/s 

105 

g-cm/sec 

Momentum 

density 

m 

Pi 

ml 
' Q * 

m 

Wt 

kg/m2-s 

10_1 

*1 

g/cm2-sec 

Permeability 

1 

henry /m 

—  x  107 

47T 

— 
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Dimensions 

Physical 

Sym- 

SI 

Conversion 

Gaussian 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

Permittivity 

e 

+2„2 
t  q 

ml 3 

1 

farad/m 

36tt  x  109 

— 

Polarization 

p 

Q 

m1'1 

coulomb/m2 

3  X  105 

statcoulomb 

/cm2 

T 

l2 

P/2t 

ml 2 

1  2 

statvolt 

Potential 

v,4> 

volt 

-  x  10  2 

Pq 

t 

3 

ml 2 

ml 2 

erg/sec 

Power 

p 

t3 

t3 

watt 

107 

Power 

density 

m 

Zt3 

m 

Ip 

watt  /  m3 

10 

erg/cm3-sec 

Pressure 

p,p 

m 

Ip 

m 

Ip 

pascal 

10 

dyne/cm2 

Reluctance 

n 

q2 

i 

ampere-turn 

4tt  x  10~9 

cm”1 

ml2 

T 

/weber 

Resistance 

R 

ml2 

tq2 

t 

i 

ohm 

1  -li 

-  x  10  11 

9 

sec/cm 

ml 3 

-  x  10-9 

9 

Resistivity 

V,P 

tq 2 

t 

ohm-m 

sec 

Thermal  con¬ 

k,  k 

ml 

ml 

watt  /  m— 

105 

erg/cm-sec- 

ductivity 

p 

t3 

deg  (K) 

deg  (K) 

Time 

t 

t 

t 

second  (s) 

1 

second  (sec) 

ml 

106 

Vector 

A 

weber/m 

gauss-cm 

potential 

tq 

l 

t 

t 

l 

t 

Velocity 

V 

m/s 

102 

cm/sec 

m 

m 

kg/ m-s 

Viscosity 

It 

It 

10 

poise 

Vorticity 

C 

1 

t 

1 

t 

s-1 

1 

sec-1 

ml 2 

ml2 

Work 

w 

t2 

P 

joule 

107 

erg 
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INTERNATIONAL  SYSTEM  (SI)  NOMENCLATURE6 


Physical 

Name 

Symbol 

Physical 

Name 

Symbol 

Quantity 

of  Unit 

for  Unit 

Quantity 

of  Unit 

for  Unit 

^length 

meter 

m 

electric 

potential 

volt 

V 

*mass 

kilogram 

kg 

electric 

ohm 

to 

*time 

second 

s 

resistance 

*current 

ampere 

A 

electric 

conductance 

siemens 

s 

*temperature 

kelvin 

K 

electric 

farad 

F 

*amount  of 
substance 

mole 

mol 

capacitance 

magnetic  flux 

weber 

Wb 

^luminous 

candela 

cd 

intensity 

magnetic 

inductance 

henry 

H 

fplane  angle 

radian 

rad 

magnetic 

tesla 

T 

fsolid  angle 

steradian 

sr 

intensity 

frequency 

hertz 

Hz 

luminous  flux 

lumen 

lm 

energy 

joule 

J 

illuminance 

lux 

lx 

force 

newton 

N 

activity  (of  a 
radioactive 

becquerel 

Bq 

pressure 

pascal 

Pa 

source) 

power 

watt 

W 

absorbed  dose 
(of  ionizing 

gray 

Gy 

electric  charge 

coulomb 

C 

radiation) 

*SI  base  unit  fSupplementary  unit 


METRIC  PREFIXES 


Multiple 

Prefix 

Symbol 

Multiple 

Prefix 

Symbol 

10_1 

deci 

d 

10 

deca 

da 

10-2 

centi 

c 

102 

hecto 

h 

CO 

1 

o 

T - 1 

milli 

103 

kilo 

k 

icr6 

micro 

106 

mega 

M 

10~9 

nano 

n 

109 

giga 

G 

10-12 

pico 

p 

1012 

tera 

T 

10~15 

femto 

f 

1015 

peta 

P 

10-18 

atto 

a 

1018 

exa 

E 

13 


PHYSICAL  CONSTANTS  (SI)7 


Physical  Quantity  Symbol  Value  Units 


Boltzmann  constant 

k 

1.3807  x  10-23 

JK"1 

Elementary  charge 

e 

1.6022  x  10~19 

C 

Electron  mass 

me 

9.1094  x  10~31 

kg 

Proton  mass 

mp 

1.6726  X  10-27 

kg 

Gravitational  constant 

G 

6.6726  X  10-11 

m3s_2kg_1 

Planck  constant 

h 

6.6261  x  10-34 

Js 

h  =  h/2it 

1.0546  x  10-34 

Js 

Speed  of  light  in  vacuum 

c 

2.9979  x  10s 

m  s_1 

Permittivity  of 

eo 

8.8542  x  10-12 

F  m-1 

free  space 

Permeability  of 

Vo 

Ait  X  10-7 

Hm'1 

free  space 

Proton/electron  mass 

mp/me 

1.8362  x  103 

ratio 

Electron  charge/mass 

e/me 

1.7588  x  1011 

Ckg-1 

ratio 

Rydberg  constant 

d  me4 

tcv  8eo  2ch3 

1.0974  x  107 

m"1 

Bohr  radius 

ao  =  eoh2  /itme2 

5.2918  x  10-11 

m 

Atomic  cross  section 

'Kao2 

8.7974  x  10-21 

m2 

Classical  electron  radius 

re  —  e2 /Ait to  me2 

2.8179  X  10-15 

m 

Thomson  cross  section 

(8it/3)re2 

6.6525  X  10-29 

m2 

Compton  wavelength  of 

h/mec 

2.4263  X  10-12 

m 

electron 

h/mec 

3.8616  x  10-13 

m 

Fine-structure  constant 

a  =  e2 /2eohc 

7.2974  x  10-3 

a-1 

137.04 

First  radiation  constant 

ci  =  2ithc2 

3.7418  x  10-16 

Wm2 

Second  radiation 

C2  =  hc/k 

1.4388  x  10-2 

mK 

constant 

Stefan-Boltzmann 

a 

5.6705  X  10"8 

W  m-2K-4 

constant 
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Physical  Quantity 

Symbol 

Value 

Units 

Wavelength  associated 

V 

O 

II 

Sr 

o 

CT> 

1.2398  X  10-6 

m 

with  1  eV 

Frequency  associated 

vo  =  e/h 

2.4180  x  1014 

Hz 

with  1  eV 

Wave  number  associated 

ko  =  e/hc 

8.0655  x  105 

m 

with  1  eV 

Energy  associated  with 

hv>o 

1.6022  x  10-19 

j 

1  eV 

Energy  associated  with 

he 

1.9864  x  10-25 

j 

1  m-1 

Energy  associated  with 

me3  /8eo2h2 

13.606 

eV 

1  Rydberg 

Energy  associated  with 

k/e 

8.6174  x  10-5 

eV 

1  Kelvin 

Temperature  associated 

e/k 

1.1604  x  104 

K 

with  1  eV 

Avogadro  number 

Na 

6.0221  x  1023 

mol-1 

Faraday  constant 

F  =  NAe 

9.6485  x  104 

C  mol-1 

Gas  constant 

R  =  NAk 

8.3145 

J  K_1mol_1 

Loschmidt’s  number 

n0 

2.6868  X  1025 

m-3 

(no.  density  at  STP) 

Atomic  mass  unit 

mu 

1.6605  x  10-27 

kg 

Standard  temperature 

To 

273.15 

K 

Atmospheric  pressure 

po  =  n0kT0 

1.0133  x  105 

Pa 

Pressure  of  1  mm  Hg 

1.3332  x  102 

Pa 

(1  torr) 

Molar  volume  at  STP 

Vo  =  RTo/po 

2.2414  X  10-2 

m3 

Molar  weight  of  air 

A/air 

2.8971  X  10~2 

kg 

calorie  (cal) 

4.1868 

J 

Gravitational 

9 

9.8067 

ms-2 

acceleration 
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PHYSICAL  CONSTANTS  (cgs)7 


Physical  Quantity 

Symbol 

Value 

Units 

Boltzmann  constant 

k 

1.3807  x  10-16 

erg/deg  (K) 

Elementary  charge 

e 

4.8032  x  10-10 

statcoulomb 

(statcoul) 

Electron  mass 

me 

9.1094  x  10~28 

g 

Proton  mass 

Trip 

1.6726  X  10-24 

g 

Gravitational  constant 

G 

6.6726  x  10"8 

dyne-cm2/g2 

Planck  constant 

h 

6.6261  X  10-27 

erg-sec 

CM 

II 

4^ 

1.0546  X  10-27 

erg-sec 

Speed  of  light  in  vacuum 

c 

2.9979  X  1010 

cm/sec 

Proton/electron  mass 

mp/me 

1.8362  x  103 

ratio 

Electron  charge/mass 

e/me 

5.2728  X  1017 

statcoul/g 

ratio 

27 r2me4 

Rydberg  constant 

T>  _ 

°°  “  ch* 

1.0974  x  105 

cm 

Bohr  radius 

ao  =  h2 /me2 

5.2918  x  10-9 

cm 

Atomic  cross  section 

7r  a02 

8.7974  x  10"17 

cm2 

Classical  electron  radius 

re  =  e2 /me2 

2.8179  X  10-13 

cm 

Thomson  cross  section 

(87r/3)re2 

6.6525  X  10-25 

cm2 

Compton  wavelength  of 

h/mec 

2.4263  X  10"10 

cm 

electron 

h/mec 

3.8616  X  10-11 

cm 

Fine-structure  constant 

a  =  e2 /he 

7.2974  x  10-3 

a-1 

137.04 

First  radiation  constant 

ci  =  27T  he2 

3.7418  x  10-5 

erg-cm2/sec 

Second  radiation 

C2  =  hc/k 

1.4388 

cm-deg  (K) 

constant 

Stefan-Boltzmann 

a 

5.6705  x  10-5 

erg/cm2- 

constant 

sec-deg4 

Wavelength  associated 

Ao 

1.2398  x  10-4 

cm 

with  1  eV 
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Physical  Quantity 

Symbol 

Value 

Units 

Frequency  associated 

^0 

2.4180  X  1014 

Hz 

with  1  eV 

Wave  number  associated 

ko 

8.0655  x  103 

cm-1 

with  1  eV 

Energy  associated  with 

1.6022  x  10-12 

erg 

1  eV 

Energy  associated  with 

1.9864  x  10-16 

erg 

1  cm-1 

Energy  associated  with 

13.606 

eV 

1  Rydberg 

Energy  associated  with 

8.6174  x  10-5 

eV 

1  deg  Kelvin 

Temperature  associated 

1.1604  x  104 

deg  (K) 

with  1  eV 

Avogadro  number 

Na 

6.0221  x  1023 

mol-1 

Faraday  constant 

F  —  N  Ae 

2.8925  X  1014 

statcoul/mol 

Gas  constant 

R  =  NAk 

8.3145  x  107 

erg/deg-mol 

Loschmidt’s  number 

n0 

2.6868  x  1019 

cm-3 

(no.  density  at  STP) 

Atomic  mass  unit 

mu 

1.6605  x  10~24 

g 

Standard  temperature 

To 

273.15 

deg  (K) 

Atmospheric  pressure 

Po  =  n0kTo 

1.0133  x  106 

dyne/cm2 

Pressure  of  1  mm  Hg 

1.3332  x  103 

dyne/cm2 

(1  torr) 

Molar  volume  at  STP 

Vo  =  RTo/po 

2.2414  X  104 

cm3 

Molar  weight  of  air 

Mair 

28.971 

g 

calorie  (cal) 

4.1868  x  107 

erg 

Gravitational 

9 

980.67 

cm/sec2 

acceleration 
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FORMULA  CONVERSION8 

Here  a  =  102  cm  m_1 ,  (3  =  107ergJ-1,  eo  =  8.8542  X  10-12Fm_1, 
H o  =  47rx  10-7  H  m_1,  c  =  (eo/^o)-1^2  =  2.9979  x  108  ms-1,  and  h  =  1.0546  x 
10-34  J  s.  To  derive  a  dimensionally  correct  SI  formula  from  one  expressed  in 
Gaussian  units,  substitute  for  each  quantity  according  to  Q  =  kQ,  where  k  is 
the  coefficient  in  the  second  column  of  the  table  corresponding  to  Q  (overbars 
denote  variables  expressed  in  Gaussian  units).  Thus,  the  formula  ao  =  h  /me 
for  the  Bohr  radius  becomes  cmo  =  (ft/3)2/[(m/3/a:2)(e2a/3/47reo)],  or  ao  = 
toh2  /  itme2 .  To  go  from  SI  to  natural  units  in  which  h  =  c  =  1  (distinguished 
by  a  circumflex),  use  Q  =  k_1Q,  where  k  is  the  coefficient  corresponding  to 
Q  in  the  third  column.  Thus  ao  =  47reofc2/[(mft/c)(e2eoftc)]  =  Ait /me2.  (In 
transforming  from  SI  units,  do  not  substitute  for  eo,  /io>  or  c.) 


Physical  Quantity 

Gaussian  Units  to  SI 

Natural  Units  to  SI 

Capacitance 

a/ Ait  eo 

eo-1 

Charge 

{ot/3/ Aiteo)1/2 

(e0hc)~1/2 

Charge  density 

(, /3/Aita5e0)1/ 2 

(e0hc)~1/2 

Current 

(a;0/47reo)1/2 

(Ho/hc)1'2 

Current  density 

(/3/Aita3e0)1/2 

(ho /hc)1/2 

Electric  field 

(A7t/3eo/a3)1/2 

(e0/hc)1/2 

Electric  potential 

(A7t/3e0/a)1/2 

(e0/hc)1/2 

Electric  conductivity 

(47T€o)-1 

eo-1 

Energy 

0 

(he)-1 

Energy  density 

/3/a3 

(he)  1 

Force 

13/a 

(he)-1 

Frequency 

1 

e-1 

Inductance 

Aiteo  /  ol 

Mo-1 

Length 

a 

1 

Magnetic  induction 

(Ait  (3/ a3  no)1/ 2 

(Hohc)~1/2 

Magnetic  intensity 

(Ait  nop/ a3)1/2 

(Ho/hc)1/2 

Mass 

(3/a2 

c/h 

Momentum 

P/a 

h-1 

Power 

(3 

(he2)-1 

Pressure 

P/a3 

(he)  1 

Resistance 

Aiteo /a 

(eo /Ho)1/2 

Time 

1 

c 

Velocity 

a 

e-1 
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MAXWELL’S  EQUATIONS 


Name  or  Description 

SI 

Gaussian 

dB 

1  dB 

Faraday’s  law 

VxE  =  — - — 
dt 

V  x  E  = - — 

c  dt 

dB 

1  dB  47t 

Ampere’s  law 

VXB=-+J 

V  x  H  =  -  — —  H  J 

c  dt  c 

Poisson  equation 

V  •  D  =  p 

V  •  D  =  4-7rp 

[Absence  of  magnetic 
monopoles] 

V  •  B  =  0 

V  •  B  =  0 

Lorentz  force  on 

q  (E  +  v  x  B) 

q  (  E  +  -v  X  B  j 

charge  q 

V  c  / 

Constitutive 

D  =  eE 

D  =  eE 

relations 

B  =  pH 

B  =  pH 

In  a  plasma,  p  «  po  =  Air  x  10  7Hm  1  (Gaussian  units:  /xml).  The 
permittivity  satisfies  e  ~  eo  =  8.8542  X  10-12Fm_1  (Gaussian:  e  ~  1) 
provided  that  all  charge  is  regarded  as  free.  Using  the  drift  approximation 
v_l  =  ExB/ B2  to  calculate  polarization  charge  density  gives  rise  to  a  dielec¬ 
tric  constant  K  =  e/eo  =  1  +  36-n-  x  109p/B2  (SI)  =  l  +  47rpc2/B2  (Gaussian), 
where  p  is  the  mass  density. 

The  electromagnetic  energy  in  volume  V  is  given  by 


W 


dV (H  B  +  E  •  D) 


dV( H  B  +  E  •  D) 


Poynting’s  theorem  is 


dW  f 

—  +  /  N  <iS  = 
J  s 


dVJ  ■  E, 


(SI) 


(Gaussian) . 


where  S  is  the  closed  surface  bounding  V  and  the  Poynting  vector  (energy  flux 
across  S )  is  given  by  N  =  E  x  H  (SI)  or  N  =  cE  x  H/47T  (Gaussian). 
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ELECTRICITY  AND  MAGNETISM 


In  the  following,  e  =  dielectric  permittivity,  /x  =  permeability  of  conduc¬ 
tor,  fi'  =  permeability  of  surrounding  medium,  a  =  conductivity,  /  =  u)/27t  = 
radiation  frequency,  Km  =  n/no  and  ne  =  e/eo-  Where  subscripts  are  used, 
‘1’  denotes  a  conducting  medium  and  ‘2’  a  propagating  (lossless  dielectric) 
medium.  All  units  are  SI  unless  otherwise  specified. 


Permittivity  of  free  space 

e0  =  8.8542  x  10-12  Fm-1 

Permeability  of  free  space 

/xo  —  47t  x  10-7Hm-1 
=  1.2566  x  10-6  Hm_1 

Resistance  of  free  space 

Ro  —  (moAo)1/2  =  376.73  Q 

Capacity  of  parallel  plates  of  area 

A,  separated  by  distance  d 

C  =  eA/d 

Capacity  of  concentric  cylinders 
of  length  Z,  radii  o,  b 

C  =  2-Kel/  In (b/a) 

Capacity  of  concentric  spheres  of 
radii  a,  b 

C  =  47t  eab/  ( b  —  a) 

Self-inductance  of  wire  of  length 

Z,  carrying  uniform  current 

L  —  jjbl 

Mutual  inductance  of  parallel  wires 
of  length  Z,  radius  a,  separated 
by  distance  d 

L  =  (pl'I/Att)  [1  +  4  ln(cZ/o)] 

Inductance  of  circular  loop  of  radius 
b,  made  of  wire  of  radius  a, 
carrying  uniform  current 

L  =  b  [ln(86/a)  -  2]  +  /i/4} 

Relaxation  time  in  a  lossy  medium 

r  —  e/cr 

Skin  depth  in  a  lossy  medium 

S  =  (2/w/xcr)1/2  =  (7r//i<7)-1/2 

Wave  impedance  in  a  lossy  medium 

Z  =  [n/(e  +  icr/a;)]1/2 

Transmission  coefficient  at 

T  =  4.22  X  10_ 4 (/ KmlKe2 /^) 1/2 

conducting  surface9 
(good  only  for  T  <§;  1) 

Field  at  distance  r  from  straight  wire 
carrying  current  /  (amperes) 

Bq  =  \il /2'kt  tesla 

=  0.21 /r  gauss  (r  in  cm) 

Field  at  distance  z  along  axis  from 

Bz=  tia2I/[2(a2  +  z2)3/2] 

circular  loop  of  radius  a 
carrying  current  I 
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ELECTROMAGNETIC  FREQUENCY/ 
WAVELENGTH  BANDS10 


Designation 

Frequency  Range 

Wavelength  Range 

Lower 

Upper 

Lower 

Upper 

ULF* 

30  Hz 

10  Mm 

VF* 

30  Hz 

300  Hz 

1  Mm 

10  Mm 

ELF 

300  Hz 

3  kHz 

100  km 

1  Mm 

VLF 

3  kHz 

30  kHz 

10  km 

100  km 

LF 

30  kHz 

300  kHz 

1  km 

10  km 

MF 

300  kHz 

3  MHz 

100  m 

1  km 

HF 

3  MHz 

30  MHz 

10m 

100  m 

VHF 

30  MHz 

300  MHz 

lm 

10  m 

UHF 

300  MHz 

3  GHz 

10  cm 

lm 

SHFf 

3  GHz 

30  GHz 

1  cm 

10  cm 

S 

2.6 

3.95 

7.6 

11.5 

G 

3.95 

5.85 

5.1 

7.6 

J 

5.3 

8.2 

3.7 

5.7 

H 

7.05 

10.0 

3.0 

4.25 

X 

8.2 

12.4 

2.4 

3.7 

M 

10.0 

15.0 

2.0 

3.0 

P 

12.4 

18.0 

1.67 

2.4 

K 

18.0 

26.5 

1.1 

1.67 

R 

26.5 

40.0 

0.75 

1.1 

EHF 

30  GHz 

300  GHz 

1  mm 

1  cm 

Submillimeter 

300  GHz 

3  THz 

100  jura 

1  mm 

Infrared 

3  THz 

430  THz 

700  nm 

100  /nm 

Visible 

430  THz 

750  THz 

400  nm 

700  nm 

Ultraviolet 

750  THz 

30  PHz 

10  nm 

400  nm 

X  Ray 
Gamma  Ray 

30  PHz 

3  EHz 

3  EHz 

100  pm 

10  nm 

100  pm 

In  spectroscopy  the  angstrom  is  sometimes  used  (lA  =  10  8  cm  =  0.1  nm). 
*The  boundary  between  ULF  and  VF  (voice  frequencies)  is  variously  defined. 
fThe  SHF  (microwave)  band  is  further  subdivided  approximately  as  shown.11 


21 


AC  CIRCUITS 


For  a  resistance  R,  inductance  L,  and  capacitance  C  in  series  with 
a  voltage  source  V  —  Vo  exp (icot)  (here  i  =  y/—l),  the  current  is  given 
by  J  =  dq/dt,  where  q  satisfies 


+  R 


dq 

dt 


=  V. 


Solutions  are  q(t)  =  qs  +  qt,  I(t )  =  Is  +  It,  where  the  steady  state  is 
Is  =  iujqs  —  VIZ  in  terms  of  the  impedance  Z  =  R  +  i(uL  —  l/u>C)  and 
It  =  dqt/dt.  For  initial  conditions  q( 0)  =  qo  =  qo  +  qs,  1(0)  =  Io>  the 
transients  can  be  of  three  types,  depending  on  A  =  R2  —  AL/C: 

(a)  Overdamped,  A  >  0 


Qt  = 

It  = 


J0  +7+90 


7+  “  7- 
7+ (A)  +  7- go) 
7+  -  7- 


exp(— 7-t) 


exp(— 7+t)  - 


—  -----  exp(-7+t), 
7+  -  7- 

7_(J0  +  7+9o) 


7+  -  7- 


exp(— 7_t), 


where  =  (R  ±  A1/2)/2L; 

(b)  Critically  damped,  A  =  0 


qt  =  [qo  +  (Jo  +  lRQo)t]  exp(-7 fit), 

It  =  [fo  -  ( Io  +  iRqohRt]  exp(-7 Rt), 


where  7 r  =  JJ/2L; 

(c)  Underdamped,  A  <  0 


Qt  = 

It  = 


7R9o  +  Jo  .  a  J  ^ 

- sm  u\t  +  qo  cos  to  it  exp(— 7nt), 

J 

r  ,  (^i2  +  7R2)<?o  +  7R^o  .  ,  ,x 

Jo  cos  u>i  t - sm(u>i  t) 

CUl 


exp(-7fit), 


Here  u>  1  =  u>o(l  —  R2C/AL)1^2 ,  where  tuo  =  ( LC)~ x^2  is  the  resonant 
frequency.  At  tu  =  ujo,  Z  =  R.  The  quality  of  the  circuit  is  Q  =  uq  L/R. 
Instability  results  when  L,  R,  C  are  not  all  of  the  same  sign. 
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DIMENSIONLESS  NUMBERS  OF  FLUID  MECHANICS 


12 


Name(s) 

Symbol 

Definition 

Significance 

Alfven, 

Karman 

Al,  Ka 

VA/V 

*  (Magnetic  force/ 
inertial  force)1/2 

Bond 

Bd 

(/ p '  - 

Gravitational  force/ 
surface  tension 

Boussinesq 

B 

V/{2  gR)1'2 

(Inertial  force/ 

gravitational  force)1/2 

Brinkman 

Br 

pV2/kAT 

Viscous  heat /conducted  heat 

Capillary 

Cp 

f*vyE 

Viscous  force/surface  tension 

Carnot 

Ca 

(T2  -T!)/t2 

Theoretical  Carnot  cycle 
efficiency 

Cauchy, 

Hooke 

Cy,  Hk 

pV2/ r  =  m2 

Inertial  force/ 

compressibility  force 

Chandra¬ 

sekhar 

Ch 

B2  L2 / pup 

Magnetic  force/dissipative 
forces 

Clausius 

Cl 

LV3  p/kAT 

Kinetic  energy  flow  rate/heat 
conduction  rate 

Cowling 

C 

(VA/V)2  =  Al2 

Magnetic  force/inertial  force 

Crispation 

Cr 

fiK/YtL 

Effect  of  diffusion/effect  of 
surface  tension 

Dean 

D 

D3/2V/u(2r)1/2 

Transverse  flow  due  to 

curvature/longitudinal  flow 

[Drag 

coefficient] 

Cd 

(p'  ~  p)Lg/ 

p'V2 

Drag  force/ inertial  force 

Eckert 

E 

V2 /cp  AT 

Kinetic  energy /change  in 
thermal  energy 

Ekman 

Ek 

(u/2QL2)1/2  = 
(Ro/Re)1/2 

(Viscous  force/Coriolis  force)1/2 

Euler 

Eu 

Ap/pV2 

Pressure  drop  due  to  friction/ 
dynamic  pressure 

Froude 

Fr 

V/(gL)1/2 

V/NL 

f  (Inertial  force/gravitational  or 
buoyancy  force)1/2 

Gay-Lussac 

Ga 

1//3AT 

Inverse  of  relative  change  in 
volume  during  heating 

Grashof 

Gr 

gL3  )3AT /u2 

Buoyancy  force/viscous  force 

[Hall 

coefficient] 

Ch 

A/ rL 

Gyrofrequency/ 
collision  frequency 

*(t)  Also  defined  as  the  inverse  (square)  of  the  quantity  shown. 
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Name(s) 

Symbol 

Definition 

Significance 

Hartmann 

H 

BL/(/xt7)1/2  = 
(Rm  Re  C)1/2 

(Magnetic  force/ 
dissipative  force) 1|/2 

Knudsen 

Kn 

X/L 

Hydrodynamic  time/ 
collision  time 

Lewis 

Le 

k/V 

^Thermal  conduction/molecular 
diffusion 

Lorentz 

Lo 

V/c 

Magnitude  of  relativistic  effects 

Lundquist 

Lu 

VoLVa/v  - 
A1  Rm 

J  x  B  force/resistive  magnetic 
diffusion  force 

Mach 

M 

V/Cs 

Magnitude  of  compressibility 
effects 

Magnetic 

Mach 

Mm 

V/VA  =  Al"1 

(Inertial  force /magnetic  force)1/2 

Magnetic 

Reynolds 

Rm 

HoLV/r) 

Flow  velocity /magnetic  diffusion 
velocity 

Newton 

Nt 

F/pL2V2 

Imposed  force/inertial  force 

Nusselt 

N 

aL/k 

Total  heat  transfer/thermal 
conduction 

Peclet 

Pe 

LV/k 

Heat  convection/heat  conduction 

Poisseuille 

Po 

D2Ap/pLV 

Pressure  force/ viscous  force 

Prandtl 

Pr 

V  /  K 

Momentum  diffusion/ 
heat  diffusion 

Rayleigh 

Ra 

gH3f3AT/vK 

Buoyancy  force/diffusion  force 

Reynolds 

Re 

LV/u 

Inertial  force/ viscous  force 

Richardson 

Ri 

(. NH/AV )2 

Buoyancy  effects/ 
vertical  shear  effects 

Rossby 

Ro 

V/2CIL  sin  A 

Inertial  force/Coriolis  force 

Schmidt 

Sc 

u/V 

Momentum  diffusion/ 
molecular  diffusion 

Stanton 

St 

a/ pepV 

Thermal  conduction  loss/ 
heat  capacity 

Stefan 

Sf 

crLT3 /k 

Radiated  heat/conducted  heat 

Stokes 

S 

v/L2f 

Viscous  damping  rate/ 
vibration  frequency 

Strouhal 

Sr 

fL/V 

Vibration  speed/flow  velocity 

Taylor 

Ta 

(2£IL2/v)2 

R1/2(AR)3/2 

■(Cl/u) 

Centrifugal  force/ viscous  force 
(Centrifugal  force/ 
viscous  force)1/2 

Thring, 

Boltzmann 

Th,  Bo 

pcp  V /  ecrT3 

Convective  heat  transport/ 
radiative  heat  transport 

Weber 

W 

pLV2/X 

Inertial  force/surface  tension 
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Nomenclature: 

B 

C  s  i  c 

Cp 

D  =  2R 
F 

f 

9 

if,  L 

k  =  pcvK 

N=(g/H)1/2 

R 

r 

rL 

T 

V 

VA  =  B/ (nop)1/2 

a 

P 

r 

AR,AV,Ap,AT 

e 

V 

k,  T> 

A 

A 

p  =  pu 

Mo 

v 

P 

p' 

S 

<7 

Cl 


Magnetic  induction 
Speeds  of  sound,  light 

Specific  heat  at  constant  pressure  (units  m2s-2  K-1) 

Pipe  diameter 

Imposed  force 

Vibration  frequency 

Gravitational  acceleration 

Vertical,  horizontal  length  scales 

Thermal  conductivity  (units  kgm-1  s~2) 

Brunt-Vaisala  frequency 

Radius  of  pipe  or  channel 

Radius  of  curvature  of  pipe  or  channel 

Larmor  radius 

Temperature 

Characteristic  flow  velocity 
Alfven  speed 

dT 

Newton’s-law  heat  coefficient,  k— —  =  a  AT 

dx 

Volumetric  expansion  coefficient,  dV/V  =  /3dT 
Bulk  modulus  (units  kgm-1  s-2) 

Imposed  differences  in  two  radii,  velocities, 
pressures,  or  temperatures 

Surface  emissivity 

Electrical  resistivity 

Thermal,  molecular  diffusivities  (units  m2s-1) 
Latitude  of  point  on  earth’s  surface 
Collisional  mean  free  path 
Viscosity 

Permeability  of  free  space 
Kinematic  viscosity  (units  m2  s-1) 

Mass  density  of  fluid  medium 

Mass  density  of  bubble,  droplet,  or  moving  object 
Surface  tension  (units  kgs-2) 

Stefan-Boltzmann  constant 
Solid-body  rotational  angular  velocity 
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SHOCKS 


At  a  shock  front  propagating  in  a  magnetized  fluid  at  an  angle  9  with 
respect  to  the  magnetic  induction  B,  the  jump  conditions  are  13,14 

(1)  pU  —  pU  =  g; 

(2)  pU2  +  p  +  Bj2 /2p  =  pU 2  +  p+  B2 / 2/x; 

(3)  pUV  -  B\\B±/p  =  pUV  -  ByBx/p; 

(4)  B||=% 

(5)  UB±  -  VB\\  =  U B±  -  V% 

(6)  \(U2  +  V2)  +  w  +  (UB2  -  V  B\\B±_)  /  ppU 

=  \  (U2  +  V2)  +  w  +  (UB2  -  VB\\  B±)/ppU. 

Here  U  and  V  are  components  of  the  fluid  velocity  normal  and  tangential  to 
the  front  in  the  shock  frame;  p  =  1/v  is  the  mass  density;  p  is  the  pressure; 
B±  =  BsinO,  J3||  =  Bcos9;  p  is  the  magnetic  permeability  (p,  =  47T  in  cgs 
units);  and  the  specific  enthalpy  is  w  =  e  +  pv,  where  the  specific  internal 
energy  e  satisfies  de  =  Tds  —  pdv  in  terms  of  the  temperature  T  and  the 
specific  entropy  s.  Quantities  in  the  region  behind  (downstream  from)  the 
front  are  distinguished  by  a  bar.  If  B  =  0,  then15 

(7)  U  -  U  =  [(p  -  p)(v  -  v)}1'2-, 

(8)  (p  -  p)(v  -  u)-1  =  g2; 

(9)  w  -  w  =  \(p  -  p)(v  +  v); 

(10)  e  -  e  =  \  (p  +  p)(v  -  v). 

In  what  follows  we  assume  that  the  fluid  is  a  perfect  gas  with  adiabatic  index 
7  =  1  +  2/n,  where  n  is  the  number  of  degrees  of  freedom.  Then  p  =  pRT/m, 
where  R  is  the  universal  gas  constant  and  m  is  the  molar  weight;  the  sound 
speed  is  given  by  Cs2  —  (dp/dp)s  =  7 pv,  and  w  =  7e  =  7P^/(7  —  1).  For  a 
general  oblique  shock  in  a  perfect  gas  the  quantity  X  =  r-1  (U /Va)2  satisfies14 

(11)  (X—  j3/a)(X— cos2  9)2  =  X  sin2  6  -£[1  +  (r  —  l)/2o:]  X  —  cos2  6  j-,  where 
r  -  p/p,  a:=i[7  +  l-(7-  l)r],  and  (3  =  Cs2/Va2  =  4tvjp/B2. 

The  density  ratio  is  bounded  by 

(12)  1  <  r  <  (7  +  1)/ (7  -  !)• 

If  the  shock  is  normal  to  B  (i.e.,  if  6  =  7r/2) ,  then 

(13)  U 2  =  (r/a)  {Cs2  +  VA2  [1  +  (1  -  j/2)(r  -  1)]}; 

(14)  U/U  =  B/B  =  r; 
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(16)  p  =  p  +  (1  -  r  1)pU2  +  (1  -  r2)B2/ 2p. 

If  0  =  0,  there  are  two  possibilities:  switch-on  shocks,  which  require  (3  <  1  and 
for  which 

(17)  U2  =  rVA2; 

(18)  U  =  Va2/U; 

(19)  B2  =  2B2(r-l){*-(3)- 

(20)  V  =  UB±/Br, 

(21)  p  =  p  +  pU2(l  -  a  +  0)(1  -  r-1), 

and  acoustic  (hydrodynamic)  shocks,  for  which 

(22)  U2  =  ( r/a)Cs 2; 

(23)  U  =  U/r; 

(24)  V  =  B±=  0; 


(25)  p  =  p  +  pU2{\  —  r-1). 

For  acoustic  shocks  the  specific  volume  and  pressure  are  related  by 

(26)  v/v  =  [(7  +  l)p  +  (7  -  1  )P]  /  [(7  -  1  )p  +  (7  +  1)p]- 
In  terms  of  the  upstream  Mach  number  M  =  U/Cs, 


(27)  p/p  =  v/v  =  U/U  =  (7  +  l)M2/[(7  -  1  )M2  +  2]; 

(28)  p/p  =  (27 M2  -  7  +  1)/ (7  +  1); 

(29)  T/T  =  [(7  -  1)M2  +  2] (2 7M2  -  7  +  l)/(7  +  1)2M2; 

(30)  M 2  =  [(7  -  1)M2  +  2]/[27M2  -  7  +  1]. 


The  entropy  change  across  the  shock  is 

(31)  A  s  =  s-s  =  cv  In  [(p/p)  (p/p)7], 

where  cv  =  R/iy  —  l)m  is  the  specific  heat  at  constant  volume;  here  R  is  the 
gas  constant.  In  the  weak-shock  limit  (M  — ►  1), 


(32)  As  — ►  Ci 


2l<y  ~  1)  (jvf2  _  i)3  ^  16t« 

3(7  +  1)  v  3(7  +  l)m 


(M  -  1) 


The  radius  at  time  t  of  a  strong  spherical  blast  wave  resulting  from  the  explo¬ 
sive  release  of  energy  £  in  a  medium  with  uniform  density  p  is 


(33)  RS  =  C0(Et2/p)1/5, 


where  Co  is  a  constant  depending  on  7.  For  7  =  7/5,  Co  =  1.033. 


FUNDAMENTAL  PLASMA  PARAMETERS 


All  quantities  are  in  Gaussian  cgs  units  except  temperature  (T,  Te,  T* ) 
expressed  in  eV  and  ion  mass  (mi)  expressed  in  units  of  the  proton  mass, 
/i  =  m,i/mp\  Z  is  charge  state;  k  is  Boltzmann’s  constant;  K  is  wavenumber; 
7  is  the  adiabatic  index;  In  A  is  the  Coulomb  logarithm. 

Frequencies 


electron  gyrofrequency 

fee  =  Uce/Z-K  =  2.80  X  106£Hz 
u>ce  =  eB/mec  =  1.76  x  107 B  rad/sec 

ion  gyrofrequency 

fci  =  u>ci/2TT  =  1.52  x  10s Z fj,-1  B  Hz 

0JC i  =  ZeB/rriiC  =  9.58  X  103Z/i-1B  rad/sec 

electron  plasma  frequency 

fpe  =  ujpe/2ir  =  8.98  X  10 3ne1^2  Hz 
upe  =  (47mee2/me)1/2 

=  5.64  x  104ne1//2  rad/sec 

ion  plasma  frequency 

fpi  =  ^L>pi/  27T 

=  2.10  x  102Z^_1/2ni1/2  Hz 

Upi  =  (47 rriiZ2e2 /mi)1/2 

=  1.32  x  103ZjU~1//27J.i1/,2rad/sec 

electron  trapping  rate 

uTe  =  ( eKE/me)1' 2 

=  7.26  x  108K'1/2E1/2  sec-1 

ion  trapping  rate 

uTi  =  ( ZeKE/mi)1/ 2 

=  1.69  x  107Z1/2AT1/2E1/:V~1/2sec_1 

electron  collision  rate 

ve  =  2.91  X  10-6ne  In  ATe-3//2  sec-1 

ion  collision  rate 

Vi  =  4.80  X  10_8Z4/i“1/2ni  In  ATi_3/2  sec-1 

Lengths 

electron  deBroglie  length 

X  =  h/(mekTe)1/2  =  2.76  x  10_8Te_1/2  cm 

classical  distance  of 

e2/kT  =  1.44  x  10-rT-1  cm 

minimum  approach 

electron  gyroradius 

T  e  =  VTe/^ce  =  2.38  Te1/2B  1  cm 

ion  gyroradius 

ri  =  VTi  / 

=  1.02  x  10 2n1/2Z~1Ti1/2B~1  cm 

plasma  skin  depth 

c/ujpe  =  5.31  X  105ne_1^2cm 

Debye  length 

\D  =  ( kT/Airne2)1' 2 

=  7.43  X  102T1/,2ti-1/2  cm 
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Velocities 

electron  thermal  velocity 
ion  thermal  velocity 
ion  sound  velocity 

Alfven  velocity 

Dimensionless 

(electron/proton  mass  ratio)1/2 

number  of  particles  in 
Debye  sphere 

Alfven  velocity /speed  of  light 

electron  plasma/gyrofrequency 
ratio 

ion  plasma/gyrofrequency  ratio 

thermal/magnetic  energy  ratio 

magnetic/ion  rest  energy  ratio 

Miscellaneous 

Bohm  diffusion  coefficient 

transverse  Spitzer  resistivity 


VTe  =  ( kTe/me)1/ 2 

=  4.19  x  107Te1//2  cm/sec 
vTi  =  ( kTi/mi)1/ 2 

=  9.79  X  10V-1/2Ti 1/2  cm/sec 
Cs  =  (7 ZkTe/mi )1/2 

=  9.79  x  lO^qZTe/p)1/2  cm/sec 
va  =  B/(47rnimi)1/2 

=  2.18  x  101  V~1/,2nj-1/213  cm/sec 

(me/mp)1/2  =  2.33  X  1CT2  =  1/42.9 
(47r/3)nAjr>3  =  1-72  x  109T3/2n_1/2 

vA/c  =  7.28  n~1/2m-1/2B 
ujpe/ojce  =  3.21  x  10-3ne1/2B-x 

ujpi/ua  =  0.137iu1/2ni1/2B-1 
/ 3  =  8-rrnkT / B2  =  4.03  X  10_11nTB“2 
B2 /8irriimiC2  =  26.5  B2 

Db  =  (ckT/16eB) 


=  6.25  X  106TB  1  cm2/sec 
r)±  =  1.15  X  10_14ZlnAT_3/2sec 
=  1.03  x  10_2Z  In  A t_3/2  ft  cm 
The  anomalous  collision  rate  due  to  low-frequency  ion-sound  turbulence  is 

v*  »  UpeW/kT  =  5.64  x  104ne1/2W' /kT sec-1, 


where  W  is  the  total  energy  of  waves  with  oj/K  <  vti- 
Magnetic  pressure  is  given  by 

Pmag  =  B2 / 8tt  =  3.98  x  106(B/Bo)2  dynes/cm2  =  3.93(B/B0)2  atm, 
where  Bo  =  10  kG  =  1  T. 

Detonation  energy  of  1  kiloton  of  high  explosive  is 

WkT  =  1012  cal  =  4.2  X  1019  erg. 
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PLASMA  DISPERSION  FUNCTION 


Definition16  (first  form  valid  only  for  ImC  >  0): 


r+°°  dt  exp  (— t2)  fl<*  .  . 

Z(()  =  n~1/2  - — — -  =  2iexp(-C2)  /  dt  exp  (-t  )  . 

«/  —  oo  «/  — oo 


Physically,  £  =  x  +  iy  is  the  ratio  of  wave  phase  velocity  to  thermal  velocity. 
Differential  equation: 

dZ  1/2  d2  Z  dZ 

—  =  -2  (1  +  C Z) ,  Z(0)  =  i7r1/2;  +  2C—  +  2Z  =  0. 

dQ  dC 

Real  argument  (y  =  0): 


Z(x)  =  exp 


Imaginary  argument  (sc  =  0): 

Z(iy)  =  i7r1/2  exp  (y2)  [1  -  erf(y)] . 


Power  series  (small  argument): 

Z( C)  =  m1/2  exp  (-C2)  -  2C  (l  -  2C2/3  +  4C4/15  -  8C6/105  +  •••). 
Asymptotic  series,  |C|  1  (Ref.  17): 

Z(0  =  iir1/2c t  exp  (-C2)  -  C_1  (l  +  1/2C2  +  3/4C4  +  15/8C6  +  •••), 
where 

[  0  y  >  |z|_1 

lyl  <  I®!-1 

I  2  y  <  -|x|_1 

Symmetry  properties  (the  asterisk  denotes  complex  conjugation): 

Z(C*)  =  -[Z(-C)]*; 


Z(C)  =  [Z( 0]  *  +  2m1/2  exp[-«*)2]  (y  >  0). 

Two-pole  approximations18  (good  for  £  in  upper  half  plane  except  when  y  < 
7r1//2a:2  exp(— x2),  x  1): 


z(  0 

Z\ 0 


0.50  +  0.81?  _  0.50  -  0.8U 
a  -  C  a*  +  C 

0.50  +  0.96i  ,  0.50  —  0.96i 

<6  -  a2  +  (>>*  +  o2 


a  =  0.51  -  0.81i; 


b  =  0.48  -  0.9H. 
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COLLISIONS  AND  TRANSPORT 


Temperatures  are  in  eV ;  the  corresponding  value  of  Boltzmann’s  constant 
is  k  =  1.60  X  10  12erg/eV;  masses  /x,  fj,'  are  in  units  of  the  proton  mass; 
ea  =  Zae  is  the  charge  of  species  a.  All  other  units  are  cgs  except  where 
noted. 

Relaxation  Rates 

Rates  are  associated  with  four  relaxation  processes  arising  from  the  in¬ 
teraction  of  test  particles  (labeled  a)  streaming  with  velocity  va  through  a 
background  of  field  particles  (labeled  (3): 


slowing  down 
transverse  diffusion 
parallel  diffusion 
energy  loss 


dv 


oc  a\f3 

~  =  Va 


dt 

d  ,  —  \2  a\/3  2 

—  (va  -  Vq)_l  =  I'j/  Va 

d  ,  —  \2  oc\/3  2 

—  (Va  -  Va)||  =  I^i  ^Va 

lLv  2  _  _U<*\0V  2 

at 


where  the  averages  are  performed  over  an  ensemble  of  test  particles  and  a 
Maxwellian  field  particle  distribution.  The  exact  formulas  may  be  written19 


v 

v 

v 

V 


<*\0 

s 

a\f3 

± 

ct\p 


<*\f3 

€ 


(1  +  ma/mp)'ilj(x0‘^l3)i'Q 

2  [(1  -  1/2 xa^)ip(xa^)  +  4>'(xa^)]  u 

['tp(xa'K(3)/xa'Kl3] 

2  [(mQ/m^)f(iQ^)  -  V,/(zQ^/3)"' 


0  ’ 


where 


ct\(3  A  2  2  \  /  2  3  2  /rki  rr> 

=  47rea  ep  Xapnp/ma  va  ;  x  =  mpva  /2kTp\ 

px 

i>(x)  =  -/=  I  dtt1/2e~*i  $'(x)  = 

^  Jo  dx 

and  Xap  =  lnAa/3  is  the  Coulomb  logarithm  (see  below).  Limiting  forms  of 
us,  u±  and  u\\  are  given  in  the  following  table.  All  the  expressions  shown 

have  units  cm3  sec-1.  Test  particle  energy  e  and  field  particle  temperature  T 
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are  both  in  eV;  fj,  =  rrii/mp  where  mp  is  the  proton  mass;  Z  is  ion  charge 
state;  in  electron-electron  and  ion-ion  encounters,  field  particle  quantities  are 
distinguished  by  a  prime.  The  two  expressions  given  below  for  each  rate  hold 

for  very  slow  ( x <C  1)  and  very  fast  ( x 1)  test  particles,  respectively. 


Electron-electron 

Jne'Ke'  ~ 
/*VAee'  « 

i/JXe//ne/Aee/  « 

Electron-ion 

ffliZ  \ei  ft 


Slow 

5.8  X  10_6T_3/2 
5.8  x  10-6T-1/2e-1 


2.9  X  10  T 


-i/2€-i 


0.23/i3/2T-3/2 


ue^i/niZ2Xei  «  2.5  x  KTV^T'  ~'~e 


\  — 4.1/2m— 1/2— 1 


_L 

/niZ2Xei  ss 


1.2  x  10-V/2T~1/2e_1- 


Ion-electron 
v?e/neZ2  Xie 
vl}e/neZ2Xie  rs  3.2  x  10 
i/^e/neZ2Aie 
Ion-ion 

,i\i' 


1.6  X  10“9Ai_1T"3/2 

-V'T-^e"1 

1.6  x  10_9/i_1T_1/2e_1 


Fast 

7.7  X  10-6e-3/2 
7.7  X  10_6e"3/2 
3.9  X  10~6Te_5/2 


3.9  x  10_6e~3/2 
7.7  X  10_6e~3/2 
2.1  x  10_9M_1Te_5/2 


1.7  X  10-V1/2e-3/2 

1.8  x  10 ~7n~1/2e~3/2 

-V/2Te"5/2 


u: 


ni,Z2Z'2Xii, 


6.8  X  10' 


-8/* 


/1/2 


1  + 


-) 

A*  ) 


-►  1.7  x  10 

-i/2 

rjp —  3/2 

-+  9.0  x  10“ 


1  1  \  M1/2 

/Li  \x' )  e3/2 


i' 


V 

_L 


n 


;/Z2Z/2Ai,/ 


1.4  x  10-V1/VlT-1/V1 


i\if 


1.8  X  10-V1/2e“3/2 


n 


;/Z2Z'2  A,,  / 


6.8  X  10-8/Li,1/2/i-1T-1/2e-1 


9.0  X  10_8/i1/2/U/_1Te  5/2 


In  the  same  limits,  the  energy  transfer  rate  follows  from  the  identity 

=  2^s  -  U±  -  I/||  , 

except  for  the  case  of  fast  electrons  or  fast  ions  scattered  by  ions,  where  the 
leading  terms  cancel.  Then  the  appropriate  forms  are 


v 


e\i 


4.2  x  10  9mZ2XPi 


e  3/2/u  1  —  8.9  X  104(/li /T)1/,2e  1  exp(  —  1836/tte/T)  sec  1 
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and 


1.8  x  10~7 Z2 Z'2 Xu/ 

^e~3^2  p1^2  /  /j,  —  1.1  (p1  /T)1^2e~1  exp(— /^e/T)j  sec  1. 


In  general,  the  energy  transfer  rate  vf^  is  positive  for  e  >  ea*  and  nega¬ 
tive  for  e  <  ea*,  where  cc*  =  (mp /ma)ea* /Tp  is  the  solution  of  if)' (x*)  = 
( ma\mp)ip(x *).  The  ratio  ea*/Tp  is  given  for  a  number  of  specific  a,  (3  in  the 
following  table: 


a\/3 

i\e  e\e,  e\p  e\D  e\T,  e\He3  e\He4 

e«* 

1.5  0.98  4.8  x  10-3  2.6  x  10~3  1.8  x  10~3  1.4  x  10-3 

When  both  species  are  near  Maxwellian,  with  T*  ^  Te,  there  are  just 
two  characteristic  collision  rates.  For  Z  =  1, 


ve  =  2.9  x  10-6nATe-3/2  sec-1 ; 

Vi  =  4.8  X  10_8nATi_3/2M_1/2sec_1. 


Temperature  Isotropization 

Isotropization  is  described  by 


dTx 

dt 


1  dTll 

2  dt 


where,  if  A  =  T±/T\\  —  1  >  0, 


2\/7t eQ  e/3  tiaAa/3  .  _2 
mQ1/2(/cT||)3/2 


-3  +  04  +  3) 


tan-^A1/2) 

AV2 


If  A  <  0,  tan  1(A1/2)/A1/2  is  replaced  by  tanh  1(— A)1^2/(— A)1/2.  For 
Tj.  «  T|,  =  T, 

z4  =  8.2  x  10_7nAT_3/2  sec-1; 

v\  =  1.9  x  10-8nAZ2/x-1/2T_3/2sec_1. 
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Thermal  Equilibration 

If  the  components  of  a  plasma  have  different  temperatures,  but  no  rela¬ 
tive  drift,  equilibration  is  described  by 

/3 


where 


p<x\P  _  18  x  1q~19  (^-a^-/3)1,/2^a2-^/32»-/3AQ/3  ^-1 


(maT/ 3  +  mpTa)3/2 
For  electrons  and  ions  with  Te  «  Ti  =  T,  this  implies 


P'V/m  =  P*Xe/ne  =  3.2  x  10_9Z2A//xT3/2cm3sec_1. 


Coulomb  Logarithm 

For  test  particles  of  mass  ma  and  charge  ea  =  Zae  scattering  off  field 
particles  of  mass  mp  and  charge  ep  =  Z^e,  the  Coulomb  logarithm  is  defined 

as  A  =  InA  =  ln^max/rmin).  Here  rm jn  is  the  larger  of  eaep/mapu2  and 
Ti/2mapu ,  averaged  over  both  particle  velocity  distributions,  where  map  = 
mamp / (mQ  +  mp)  and  u  =  va  —  vp ;  rmax  =  (47t  n7e72/fcT7)_1//2,  where 

the  summation  extends  over  all  species  7  for  which  u2  <  vt^2  =  fcT7/m7.  If 
this  inequality  cannot  be  satisfied,  or  if  either  uu)Ca.  <  rmax  or  uucp  < 
rmaxi  the  theory  breaks  down.  Typically  A  «  10-20.  Corrections  to  the  trans¬ 
port  coefficients  are  0(A-1);  hence  the  theory  is  good  only  to  ~  10%  and  fails 
when  A  rsj  1  # 

The  following  cases  are  of  particular  interest: 

(a)  Thermal  electron-electron  collisions 

Aee  =23 -In  (ne1/2Te_3/2), 

=  24 -In  (ne1/2Te_1), 

(b)  Electron-ion  collisions 

Aei  =  Aie  =  23  -  In  (ne1/2ZTe-3/2)  , 

=  24 -In  (ne17^;1), 

=  30  -  In  (ni1/2Ti~3/2Z2n~1), 


Te  <  10  eV; 

Te  >  10  eV. 

Time/rrii  <  Te  <  10Z2eV; 
Time/rrii  <  10 Z2  eV  <  Te 
Te  <  TiZme/rrii. 
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(c)  Mixed  ion-ion  collisions 


zz 


'Jj^  +  T)  /n,Z2  ryZ/2\ 

i'  +  fl'Ti  \  Ti  Ttr  ) 


'2\  1/2' 


(d)  Counterstreaming  ions  (relative  velocity  vd  =  Pd c)  in  the  presence  of 
warm  electrons,  kTi/rrii^kT^  /m^  <  uo2  <  kTe/me 


ZZ\^±T)  /ne\1/2 

W'Pd2  \Te 


Fokker-Planck  Equation 

D/a  _  a/Q 


Dt  dt 


+  v  •  V/a  +  F  •  Vv/ 


«  =  far 

V  dt 


coll 


where  F  is  an  external  force  field.  The  general  form  of  the  collision  integral  is 
(dfa/dt) con  =  -  Vv  •  JqV/3,  with 


JqX/3  =  2ttA 


2  2  C 

eot  e{3  I  ,3  t  ,  2  /  \  —3 

a/3 -  /  d  V  [U  I  —  UU)U 

ma  J 

■  | — r  <v)vv,/V)  -  -/(»')Vv/"(v)| 

^  m/3  ma  J 

(Landau  form)  where  u  =  v'  —  v  and  /  is  the  unit  dyad,  or  alternatively, 

|/“ (v)V„H(v)  -  ivv  •  [/“  (v)VvV„G(v)]  |  , 


J^‘s=4^aS^f^ 

ma2 


where  the  Rosenbluth  potentials  are 


■/ 


G(v)  =  /  f^(v')ud3v' 


#  (v)  =  1  + 


mr 


m/3 


1d3v'. 
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If  species  a  is  a  weak  beam  (number  and  energy  density  small  compared  with 
background)  streaming  through  a  Maxwellian  plasma,  then 


ja\/3  = - 

ma  +  mp  3  2  II 

-  j-'f'5  (^2/-w  )v,r. 


B-G-K  Collision  Operator 

For  distribution  functions  with  no  large  gradients  in  velocity  space,  the 
Fokker-Planck  collision  terms  can  be  approximated  according  to 

=  Vee(Fe  -  fe )  +  Vei(Fe  -  fe)\ 


=  yie{Fi  -  fi)  +  l>h (Fi  -  fi). 

The  respective  slowing-down  rates  given  in  the  Relaxation  Rate  section 

above  can  be  used  for  ya0i  assuming  slow  ions  and  fast  electrons,  with  e  re¬ 
placed  by  Ta.  (For  vee  and  ua,  one  can  equally  well  use  v± ,  and  the  result 
is  insensitive  to  whether  the  slow-  or  fast-test-particle  limit  is  employed.)  The 
Maxwellians  Fa  and  Fa  are  given  by 


F(x  ^a 


Fa  —  TIq, 


/  ma  \ 

\  2nkTa  ) 


3/2 


exp 


}• 


where  na ,  va  and  Ta  are  the  number  density,  mean  drift  velocity,  and  effective 
temperature  obtained  by  taking  moments  of  fa  •  Some  latitude  in  the  definition 
of  Ta  and  vQ  is  possible;20  one  choice  is  Te  =  Tj,  T*  =  Te,  ve  =  v*,  v*  =  ve. 

Transport  Coefficients 

Transport  equations  for  a  multispecies  plasma: 


da  na 
dt 


+  7la  V 


Va  =  0; 


Tl  Q- 


dQVq 
dt 


—  Vpa  ^  *  F i a  Za^Ha 


E  H — va  x  B 

c 


+  R-a! 
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3  dakT 

-ria — ~7r~  +PaV  •  va  =  -V  •  qa  -  Pa  :  Vva  +  Qa. 

2  at 

Here  da / dt  =  d/dt  +  va  -V;  pa  —  nakTa,  where  k  is  Boltzmann’s  constant; 
Rq  =  ^Rtt)3  and  Q&  =  '^//3Qa0,  where  Ra/3  and  Qap  are  respectively 

the  momentum  and  energy  gained  by  the  ath  species  through  collisions  with 
the  /3th;  Pa  is  the  stress  tensor;  and  qQ  is  the  heat  flow. 

The  transport  coefficients  in  a  simple  two-component  plasma  (electrons 
and  singly  charged  ions)  are  tabulated  below.  Here  ||  and  _L  refer  to  the  di¬ 
rection  of  the  magnetic  field  B  =  b B\  u  =  ve  —  is  the  relative  streaming 
velocity;  ne  =  rii  =  n;  j  =  -neu  is  the  current;  uce  =  1.76  x  107Bsec-1  and 
wCi  =  ( me/rrii)u!ce  are  the  electron  and  ion  gyrofrequencies,  respectively;  and 
the  basic  collisional  times  are  taken  to  be 


re 


3  y^(fcTe)3/2 
4\/2tt  nXe4 


rp  3/2 

=  3.44  x  105-^—  sec, 
nX 


where  A  is  the  Coulomb  logarithm,  and 


r% 


3v/m7(fcTi)3^2 

4\flxn  Ae4 


=  2.09  x  107 


rp  3/2 

n\ 


sec. 


In  the  limit  of  large  fields  (u;caTa  1,  a  =  i,  e)  the  transport  processes  may 
be  summarized  as  follows:21 


momentum  transfer 
frictional  force 

electrical 

conductivities 


Hei  —  Hie  —  H  —  —f“  H.71  ? 

Ru  =  ne(j||/<7||  +  jx/o-J.); 

2 

c ||  =  1.96crj_;  cr±  =  ne  re/me\ 


thermal  force 

ion  heating 

electron  heating 
ion  heat  flux 

ion  thermal 
conductivities 
electron  heat  flux 

frictional  heat  flux 


R t=  —0.71nV|i  (kTe)  — 


3  n 


7”e 


b  X  Vj, (fcTe); 


Qi 


3 me  nk 


(Tc-T<); 


mi  T  e 

Qe  =  “ Qi  -  R  •  u; 

Qz  =  ftn  V ||  {kTi )  ftj_  V _l  (kTi)  -j-  ft/\ b  X  V  j_ (kTi)] 


i  0  ^nkTin 
ft  11  =  3.9 - 

11  mi 

qe  =qu +  9Lt; 


k±  = 


2nkTi 
miojjri  ’ 


5  nkTi 


ka  —  - ; 

2  mnvci 


0.71nfcTeui|  + 


3  nkTe 

2  Wce^e 


b  X  u_l; 
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thermal  gradient 
heat  flux 

electron  thermal 
conductivities 


qr  =  -«jjV||  (kTe)  -  Ke±V ±(kTe)  -  KeA b  x  Vx(fcTe); 


ac  ii  —  3.2 


nkTere 


mf 


;  /cl  =4.7 


nkTp 


meUJ^Te 


5 


5  nkTe 
2roew0 


stress  tensor  (either  PXx  =  —~~(Wxx  +  Wyy )  —  -^(Wxx  -  W^y)  -  r)3Wxy-, 
species)  2  2. 

Pyy=  ~^(WXX  +  »V„)  +  y(W**  -  "W  +18W',,; 

773 

Pxy—  Pyx  =  —f)\WXy  ~\  ~ (VFxx  —  Wyy)\ 

Pxz  —  Pzx  =  V%  Wxz  t74Wy*; 

Pyz~~  Pzy  =  'H2^^yz  "I-  ^74W^xz! 

Pzz  =  VoW^zz 

(here  the  z  axis  is  defined  parallel  to  B); 

3  nkTi  i  6nkTi 

~2~~  >  *?2  = 


ion  viscosity 

vl 

=  0.96nfcTjTi; 

v\  = 

% 

nkTi  i 

nkTi 

v3 

“  2Wc(  '  *  " 

W ci 

electron  viscosity 

Vo 

=  0.73  nkTeTe‘, 

vt  = 

5w  ?Ti  ’ 


nkTe  e  o  n 
;  Vi  =  2-0 


nkTe 

O  > 


*?3 


nkTe  e  nkTe 
V4  = 


2cUCe  Wee 

For  both  species  the  rate-of-strain  tensor  is  defined  as 


Wjk  = 


dv 


dvk 


H - — - <57-feV  -v. 

chc/e  3 


When  B  =  0  the  following  simplifications  occur: 

Ru  =  nej/<7ii;  Rt  =  -0.71nV(fcTe);  q7  = -/CnV(fcTi); 


q®  =  0.71nfcTeu;  q^  = -K^V(fcTe);  Pjfc  =  -VoWjk. 

For  u>ceTe  1  ojCiT i ,  the  electrons  obey  the  high-field  expressions  and  the 

ions  obey  the  zero-field  expressions. 

Collisional  transport  theory  is  applicable  when  (1)  macroscopic  time  rates 
of  change  satisfy  d/dt  1  /r,  where  r  is  the  longest  collisional  time  scale,  and 
(in  the  absence  of  a  magnetic  field)  (2)  macroscopic  length  scales  L  satisfy  L 
l,  where  l  =  vt  is  the  mean  free  path.  In  a  strong  field,  wcer  1,  condition 
(2)  is  replaced  by  Ly  l  and  L±  y/lrl  ( L±  >  re  in  a  uniform  field), 
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where  L y  is  a  macroscopic  scale  parallel  to  the  field  B  and  Lj_  is  the  smaller 
of  B/\V±B\  and  the  transverse  plasma  dimension.  In  addition,  the  standard 
transport  coefficients  are  valid  only  when  (3)  the  Coulomb  logarithm  satisfies 
A  1;  (4)  the  electron  gyroradius  satisfies  re  ^  Xp,  or  87memec2  >>  J32;  (5) 
relative  drifts  u  =  va  —  v/3  between  two  species  are  small  compared  with  the 
thermal  velocities,  i.e.,  u2  <C  kTa/ma,  kTp/mp]  and  (6)  anomalous  transport 
processes  owing  to  microinstabilities  are  negligible. 

Weakly  Ionized  Plasmas 

Collision  frequency  for  scattering  of  charged  particles  of  species  a  by 
neutrals  is 

Va.  =  nocrf^0  (kTa  /rria)1  /2 , 

where  no  is  the  neutral  density  and  crf^0  is  the  cross  section,  typically  ~ 

5  X  10~15  cm2  and  weakly  dependent  on  temperature. 

When  the  system  is  small  compared  with  a  Debye  length,  L  <  Ap,  the 
charged  particle  diffusion  coefficients  are 

Da  =  kTa/ma  Va, 

In  the  opposite  limit,  both  species  diffuse  at  the  ambipolar  rate 

HiDe  —  neDi  ( Ti  +  Te)DiDe 

Da  —  -  =  - , 

Mi  —  Me  TiDe  +  TeDi 


where  =  ea/mava  is  the  mobility.  The  conductivity  cra  satisfies  aa  = 
f^a  &a  Ma  • 

In  the  presence  of  a  magnetic  field  B  the  scalars  m  and  a  become  tensors, 


JQ  =  <ra  E  =  <rjj* E||  +  cr“  E±  +  <t“E  X  b, 


where  b  =  B /B  and 


a  2  / 

& ||  —  YI'Ol&ol  f  TTtctVoc] 

ot  oc  2//  2.  2\ 

ax  =  a ||  Va  /{Va  +  WCQ); 

O A  =  cr ||  VaOJca/{va  +  ^ca). 

Here  <t_l  and  crA  are  the  Pedersen  and  Hall  conductivities,  respectively. 
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APPROXIMATE  MAGNITUDES 
IN  SOME  TYPICAL  PLASMAS 


Plasma  Type 

n  cm  3 

T  eV 

OJpe  sec  1 

A  d  cm 

00 

Q 

£ 

Uei  sec  1 

Interstellar  gas 

1 

1 

6  x  104 

7  x  102 

4  X  10s 

7  x  10-5 

Gaseous  nebula 

h- 1 

o 

CO 

1 

2  x  106 

20 

8  x  106 

6  x  10-2 

Solar  Corona 

109 

102 

2  x  109 

2  x  10-1 

8  x  106 

60 

Diffuse  hot  plasma 

1012 

102 

6  x  1010 

7  x  10-3 

4  x  105 

40 

Solar  atmosphere, 
gas  discharge 

1014 

1 

6  x  1011 

7  x  10-5 

40 

2  x  109 

Warm  plasma 

1014 

10 

6  x  1011 

2  x  10~4 

8  x  102 

107 

Hot  plasma 

1014 

102 

6  x  1011 

7  x  10-4 

4  x  104 

4  x  106 

Thermonuclear 

plasma 

1015 

104 

2  x  1012 

2  x  10-3 

8  x  106 

5  x  104 

Theta  pinch 

1016 

102 

6  x  1012 

7  x  10-5 

4  x  103 

3  x  10s 

Dense  hot  plasma 

h-1 

o 

00 

102 

6  x  1013 

7  x  10-6 

4  x  102 

2  x  1010 

Laser  Plasma 

102° 

102 

6  x  1014 

7  x  10-7 

40 

2  x  1012 

The  diagram  (facing)  gives  comparable  information  in  graphical  form. 
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IONOSPHERIC  PARAMETERS23 


The  following  tables  give  average  nighttime  values.  Where  two  numbers 
are  entered,  the  first  refers  to  the  lower  and  the  second  to  the  upper  portion 
of  the  layer. 


Quantity 

E  Region 

F  Region 

Altitude  (km) 

90-160 

160-500 

Number  density  (m-3) 

1.5  x  10lo-3.0  X  1010 

5  X  1010-2  x  1011 

Height-integrated  number 
density  (m-2) 

9  X  1014 

4.5  X  1015 

Ion-neutral  collision 
frequency  (sec-1) 

2  x  103-102 

0.5-0.05 

Ion  gyro-/collision 
frequency  ratio 

0.09-2.0 

4.6  x  102-5.0  x  103 

Ion  Pederson  factor 

«i/(  1  +  «a2) 

0.09-0.5 

2.2  x  10-3-2  x  10-4 

Ion  Hall  factor 

Ki2 /{l  +  Ki2) 

8  x  10-4-0.8 

1.0 

Electron-neutral  collision 
frequency 

1.5  x  104-9.0  x  102 

80-10 

Electron  gyro-/collision 
frequency  ratio  Ke 

4.1  x  102-6.9  x  103 

7.8  X  104-6.2  x  105 

Electron  Pedersen  factor 

Ke/(  1  +  Ke2) 

2.7  X  10-3-1.5  X  10-4 

10-5-1.5  x  10-6 

Electron  Hall  factor 
/Ce2/(1  +  «e2) 

1.0 

1.0 

Mean  molecular  weight 

28-26 

22-16 

Ion  gyrofrequency  (sec-1) 

180-190 

230-300 

Neutral  diffusion 

coefficient  (m2sec-1) 

30-5  x  103 

105 

The  terrestrial  magnetic  field  in  the  lower  ionosphere  at  equatorial  lati¬ 
tudes  is  approximately  Bo  =  0.35  XlO-4  tesla.  The  earth’s  radius  is  Re  =  6371 
km. 
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SOLAR  PHYSICS  PARAMETERS24 


Parameter 

Symbol 

Value 

Units 

Total  mass 

Mq 

1.99  x  1033 

g 

Radius 

Rq 

6.96  x  1010 

cm 

Surface  gravity 

90 

2.74  x  104 

cms-2 

Escape  speed 

Voo 

6.18  X  107 

cm  s-1 

Upward  mass  flux  in  spicules 

— 

1.6  x  10-9 

g cm”2  s-1 

Vertically  integrated  atmospheric  density 

— 

4.28 

gem-2 

Sunspot  magnetic  field  strength 

-Bmax 

2500-3500 

G 

Surface  effective  temperature 

To 

5770 

K 

Radiant  power 

T0 

3.83  X  1033 

ergs-1 

Radiant  flux  density 

T 

6.28  x  1010 

ergcm-2s-1 

Optical  depth  at  500  nm,  measured 

T5 

0.99 

— 

from  photosphere 

Astronomical  unit  (radius  of  earth’s  orbit) 

AU 

1.50  x  1013 

cm 

Solar  constant  (intensity  at  1  AU) 

/ 

1.36  x  106 

erg  cm-2  s-1 

Chromosphere  and  Corona25 


Parameter  (Units) 

Quiet 

Sun 

Coronal 

Hole 

Active 

Region 

Chromospheric  radiation  losses 
(erg  cm-2  s-1) 

Low  chromosphere 

2  x  106 

2  X  106 

£  107 

Middle  chromosphere 

2  X  106 

2  x  106 

107 

Upper  chromosphere 

3  X  105 

3  x  105 

2  x  106 

Total 

4  x  106 

4  x  106 

>  2  x  107 

Transition  layer  pressure  (dyne cm-2) 

0.2 

0.07 

2 

Coronal  temperature  (K)  at  1.1  R© 

1. 1-1.6  x  106 

106 

2.5  x  106 

Coronal  energy  losses  (erg cm-2  s-1) 

Conduction 

2  x  105 

6  x  104 

105-107 

Radiation 

105 

104 

5  x  106 

Solar  Wind 

<  5  x  104 

7  x  105 

<  105 

Total 

3  x  105 

8  x  105 

107 

Solar  wind  mass  loss  (gem-2  s-1) 

<  2  x  10-11 

2  x  10-1° 

<  4  x  10-11 
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THERMONUCLEAR  FUSION 


26 


Natural  abundance  of  isotopes: 

hydrogen  nD/nH  =  1-5  X  10-4 

helium  nHe3/nHe4  =  x  10-6 

lithium  nLi6/nLi7  =  0-08 


Mass  ratios: 


me/mD  =  2.72  x  10-4  =  1/3670 
(me/m£>)1//2=  1.65  X  10-2  =  1/60.6 
me/mT  =  1.82  x  10-4  =  1/5496 
(rrie/mT)1^2  =  1-35  X  10-2  =  1/74.1 

Absorbed  radiation  dose  is  measured  in  rads:  1  rad  =  102ergg-1.  The  curie 
(abbreviated  Ci)  is  a  measure  of  radioactivity:  1  curie  =  3.7xl010  counts  sec-1. 

Fusion  reactions  (branching  ratios  are  correct  for  energies  near  the  cross  section 
peaks;  a  negative  yield  means  the  reaction  is  endothermic):27 

T(l.OlMeV)  +  p(3.02  MeV) 


(la) 

D 

+ 

D 

(lb) 

(2) 

D 

+ 

T 

(3) 

D 

+ 

He3 

(4) 

T 

+ 

T 

(5a) 

He3 

+  T 

(5b) 

(5c) 

(6) 

P 

+ 

Li6 

(7a) 

P 

+ 

Li7 

(7b) 

(8) 

D 

+ 

Li6 

(9) 

P 

+ 

B11 

(10) 

n 

+ 

Li6 

50% 

50% 


51% 


43% 


6% 


►He3  (0.82  MeV)  +  n(2.45MeV) 
-►He4  (3.5  MeV)  +  n(14.1MeV) 
-►He4  (3.6  MeV)  +  p(14.7MeV) 
-►He4  +  2n  +  11.3  MeV 
-►He4  +  p  +  n  +  12.1  MeV 
-►He4  (4.8  MeV)  +  D(9.5MeV) 
-►He5  (2.4  MeV)  +  p(11.9MeV) 


20% 


80% 


►  He4(1.7  MeV)  +  He3(2.3MeV) 

►  2  He4  +  17.3  MeV 

►  Be7  +  n  —  1.6  MeV 
-►2He4  +  22.4  MeV 
-►3  He4  +  8.7  MeV 

-►He4  (2.1  MeV)  +  T(2.7MeV) 


The  total  cross  section  in  barns  (1  barn  =  10-24cm2)  as  a  function  of  E,  the 
energy  in  keV  of  the  incident  particle  [the  first  ion  on  the  left  side  of  Eqs. 
(l)-(5)],  assuming  the  target  ion  at  rest,  can  be  fitted  by28 

A5  +  [(A4  -  A3E)2  +  l]_1  A2 


<7 t  (F1) 


E  [exp (AiE-1/2)  -  1 
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where  the  Duane  coefficients  Aj  for  the  principle  fusion  reactions  are  as  follows: 


1 

D-D 

(la) 

D-D 

(lb) 

D-T 

(2) 

D-He3 

(3) 

T-T 

(4) 

T-He3 

(5a-c) 

Ax 

47.88 

45.95 

89.27 

38.39 

123.1 

^2 

372 

482 

50200 

25900 

448 

11250 

^3 

1.368  x  10-2 

3.98  X  10-3 

1.02  x  10-3 

0 

a4 

1.177 

1.076 

1.297 

2.09 

0 

A5 

409 

647 

0 

0 

Reaction  rates  av  (in  cm3  sec  1),  averaged  over  Maxwellian  distributions: 


Temperature 

(keV) 

D-D 

(la  -|-  lb) 

D-T 

(2) 

D-He3 

(3) 

T-T 

(4) 

T-He3 

(5a-c) 

1.0 

1.5  x  10-22 

5.5  x  10-21 

10-26 

3.3  x  10-22 

10_ 28 

2.0 

5.4  x  10-21 

2.6  x  10~19 

1.4  x  10-23 

7.1  x  10-21 

10-25 

5.0 

1.8  x  10~19 

1.3  X  10-17 

6.7  X  10"21 

1.4  x  10-19 

2.1  x  10-22 

10.0 

1.2  x  10-18 

1.1  x  10~16 

2.3  x  10-19 

7.2  x  10-19 

1.2  x  10_2° 

20.0 

5.2  x  10-18 

4.2  x  10-16 

3.8  x  10-18 

2.5  x  10~18 

2.6  X  10-19 

50.0 

2.1  x  10-17 

8.7  x  10-16 

5.4  x  10-17 

8.7  x  10-18 

5.3  X  10-18 

100.0 

4.5  X  10-17 

8.5  x  10-16 

1.6  x  10-16 

1.9  x  10-17 

2.7  X  10-17 

200.0 

8.8  X  10-17 

6.3  x  10-16 

2.4  x  10-16 

4.2  x  10-17 

9.2  x  10-17 

500.0 

1.8  X  10-16 

3.7  x  10-16 

2.3  x  10~16 

8.4  x  10-17 

2.9  X  10-16 

1000.0 

2.2  x  10-16 

2.7  x  10-16 

1.8  x  10-16 

8.0  x  10-17 

5.2  x  10~16 

For  low  energies  (T  <  25  keV)  the  data  may  be  represented  by 

(av)DD  =  2.33  X  10_14T_2/3exp(-18.76T_1/3)cm3sec_1; 


Hdt  =  3.68  X  10“12T_2/3exp(-19.94T_1/3)cm3sec  \ 
where  T  is  measured  in  keV. 

The  power  density  released  in  the  form  of  charged  particles  is 

Pdd  —  3.3  X  lO~lznD2 (ov)dd  watt  cm-3  (including  the  subsequent 

D-T  reaction); 

Pdt  —  5.6  x  10  13nD«'T(<7f)jr)T  watt  cm  3; 

^DHe3  =  2-9  X  10  12nj->nHe3 (o‘u)£)He3  wattcm 
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RELATIVISTIC  ELECTRON  BEAMS 


Here  7  =  (1  —  /32)-1/2  is  the  relativistic  scaling  factor;  quantities  in 
analytic  formulas  are  expressed  in  SI  or  cgs  units,  as  indicated;  in  numerical 
formulas,  I  is  in  amperes  (A),  B  is  in  gauss  (G),  electron  linear  density  N  is 
in  cm-1,  and  temperature,  voltage  and  energy  are  in  MeV;  =  vz/c\  k  is 
Boltzmann’s  constant. 

Relativistic  electron  gyroradius: 

re  =  ^-(72  -  l)1/2  (cgs)  =  1.70  X  103(72  -  1)1/2£_1  cm. 
eB 

Relativistic  electron  energy: 


W  =  me2 7  =  O.5II7  MeV. 


Bennett  pinch  condition: 

I2  =  2 Nk(Te  +  Ti)c 2  (cgs)  =  3.20  x  10_4iV(Te  +  T<)  A2. 
Alfven-Lawson  limit: 

I A  =  (mc3/e)/3z 7  (cgs)  =  (4-7rmc/ n0e)/3z-y  (SI)  =  1.70  X  104/3Z7  A. 

The  ratio  of  net  current  to  I  a  is 

I  _  v 
I A  7  ’ 

Here  v  =  Nre  is  the  Budker  number,  where  re  =  e2/mc2  =  2.82  X  10-13cm 
is  the  classical  electron  radius.  Beam  electron  number  density  is 

nb  =  2.08  X  10 8J/3~1  cm-3, 

where  J  is  the  current  density  in  A  cm-2.  For  a  uniform  beam  of  radius  a  (in 
cm), 

nb  =  6.63  x  10 7 Ia~2 f3_1  cm-3, 

and 

2 re  __  z/ 
a  7 
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Child’s  law:  (non-relativistic)  space-charge-limited  current  density  between 
parallel  plates  with  voltage  drop  V  (in  MV)  and  separation  d  (in  cm)  is 


J  =  2.34  x  103V3/2d~2  A  cm-2. 


The  saturated  parapotential  current  (magnetically  self-limited  flow  along  equi- 
potentials  in  pinched  diodes  and  transmission  lines)  is29 

Ip  =  8.5  x  103G7  In  [7  +  (7s  -  1)1/2]  A, 

where  G  is  a  geometrical  factor  depending  on  the  diode  structure: 
w  for  parallel  plane  cathode  and  anode 

of  width  w,  separation  d; 

-1 

for  cylinders  of  radii  R\  (inner)  and  R2  (outer); 

for  conical  cathode  of  radius  Rc ,  maximum 
do  separation  do  (at  r  =  Rc )  from  plane  anode. 

For  (3  —*  0  (7  — ►  1),  both  I  a  and  Ip  vanish. 

The  condition  for  a  longitudinal  magnetic  field  Bz  to  suppress  filamentation 
in  a  beam  of  current  density  J  (in  A  cm-2)  is 

Bz  >47/3z(7J)1/2  G. 

Voltage  registered  by  Rogowski  coil  of  minor  cross-sectional  area  A ,  n  turns, 
major  radius  a,  inductance  L,  external  resistance  R  and  capacitance  C  (all  in 
SI): 


vjt  — 


27rd 


G  = 


In* 

Ri 


C  _ 


externally  integrated  V  =  (1/ RC)(nA/j,oI /27va)', 

self-integrating  V  =  (R/ L)(nA/j,oI /2ira)  =  Rl/n. 

X-ray  production,  target  with  average  atomic  number  Z  (V  ^  5  MeV): 

77  =  x-ray  power/beam  power  =  7  X  10  4ZV. 

X-ray  dose  at  1  meter  generated  by  an  e-beam  depositing  total  charge  Q 
coulombs  while  V  >  0.84Vmax  in  material  with  charge  state  Z: 

D  =  150Vmfx8QZ1/2  rads. 
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BEAM  INSTABILITIES 


30 


Name 

Conditions 

Saturation  Mechanism 

Electron- 

electron 

Vd  ^  Vcjf  5  3  —  1)  2 

Electron  trapping  until 

Vej  ~  Vd 

Buneman 

Vd  >  ( M/m)1/3Vi , 
Vd>Ve 

Electron  trapping  until 

Ve  ~  Vd 

Beam-plasma 

Vb  >  (np/nb)1/3Vb 

Trapping  of  beam  electrons 

Weak  beam- 
plasma 

Vb  <  ( np/nb)1/3Vb 

Quasilinear  or  nonlinear 
(mode  coupling) 

Beam-plasma 

(hot-electron) 

Ve>Vb>  Vb 

Quasilinear  or  nonlinear 

Ion  acoustic 

Te  >  Ti,  Vd  »  Cs 

Quasilinear,  ion  tail  form¬ 
ation,  nonlinear  scattering, 
or  resonance  broadening. 

Anisotropic 

temperature 

(hydro) 

Te±  >  2Te|| 

Isotropization 

Ion  cyclotron 

Vd  >  20 Vi  (for 

Te  &  Ti) 

Ion  heating 

Beam-cyclotron 

(hydro) 

vd  >  cs 

Resonance  broadening 

Modified  two- 
stream  (hydro) 

Vd  <  (1  +  /3)1/2Va, 

vd  >  cs 

Trapping 

Ion-ion  (equal 
beams) 

U  <  2(1  +  /3)1/2Va 

Ion  trapping 

Ion-ion  (equal 
beams) 

U  <  2  Cs 

Ion  trapping 

For  nomenclature,  see  p.  50. 
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Name 

Parameters  of  Most  Unstable  Mode 

Wave 

Group 

Growth  Rate 

Frequency 

Number 

Velocity 

Electron- 

1 

-We 

0 

ve 

0.9  — 

0 

electron 

2 

vd 

Buneman 

/my/3 

0.7  —  uje 

n  ,  (m\^z 

0.4  (  —  )  u>e 

Ve 

\Vd 

\m) 

\m) 

Vd 

3 

Beam-plasma 

o.t  yy/3., 

\nj 

ve  — 

0Je 

Vb 

1* 

/niX1'3 

0.4  I  I  Ve 

\np  J 

Weak  beam- 

ve 

3Ve2 

plasma 

2 np  \VbJ 

Vb 

vb 

Beam-plasma 

(hot-electron) 

(nb\1/2  Ve 

Up/  VbUe 

vb 

-=-U>e 

Ve 

XD 

vb 

Ion  acoustic 

[m) 

Vi 

Xd' 

Cs 

Anisotropic 

fie 

Ve  COS  6  rsj  Qe 

rl1 

Ve± 

temperature 

(hydro) 

Ion  cyclotron 

O.l^i 

1.2  ft* 

rt_1 

-3* 

Beam-cyclotron 

0.7  Cle 

nfle 

0.7\d 

Z  vd; 

(hydro) 

£  cs 

Modified  two- 

1 

— 

0.9  £2h 

1.7—7- 

\vd 

stream 

(hydro) 

2 

vd 

YIh 

1-2-rf 

2 

Ion-ion  (equal 

0.4  £lfj 

0 

0 

beams) 

U 

Ion-ion  (equal 

0Au>i 

0 

Vi 

1.2-77 

0 

beams) 

u 

For  nomenclature,  see  p.  50. 
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In  the  preceding  tables,  subscripts  e,  i,  d,  b ,  p  stand  for  “electron,”  “ion,” 
“drift,”  “beam,”  and  “plasma,”  respectively.  Thermal  velocities  are  denoted 
by  a  bar.  In  addition,  the  following  are  used: 


m 

electron  mass 

re,ri 

gyroradius 

M 

ion  mass 

plasma/magnetic  energy 

V 

velocity 

density  ratio 

T 

temperature 

VA 

Alfven  speed 

Tie  >  Tl% 

number  density 

gyrofrequency 

n 

harmonic  number 

Cin 

hybrid  gyrofrequency, 

Cs  =  ( Te/M )1/2 

ion  sound  speed 

Off2  =  fieOi 

plasma  frequency 

U 

relative  drift  velocity  of 

A  D 

Debye  length 

two  ion  species 

LASERS 


System  Parameters 

Efficiencies  and  power  levels  are  approximately  state-of-the-art  (1990). 31 


Type 

Wavelength 

Efficiency 

Power  levels  available  (W) 

(pm) 

Pulsed 

CW 

co2 

10.6 

0.01-0.02 

(pulsed) 

>  2  x  1013 

>  105 

CO 

5 

0.4 

>  109 

>  100 

Holmium 

2.06 

0.03f-0.1t 

>  107 

30 

Iodine 

1.315 

0.003 

>  1012 

- 

Nd-glass, 

YAG 

1.06 

0. 001-0. 06f 
>  0.1$ 

~  1014  (ten- 
beam  system) 

1-103 

*Color  center 

1-4 

10~3 

>  106 

1 

*Vibronic  (Ti 
Sapphire) 

0.7-0. 9 

>  0.1  X  rjp 

106 

1-5 

Ruby 

0.6943 

<  10-3 

1010 

1 

He-Ne 

0.6328 

h-» 

O 

1 

4^ 

- 

1-50X10-3 

*Argon  ion 

0.45-0.60 

10-3 

5  x  104 

1-20 

*OPO 

0.4-9. 0 

>  0.1  x  r]p 

o 

rH 

1-5 

n2 

0.3371 

0.001-0.05 

105-106 

- 

*Dye 

0.3-1. 1 

10-3 

>  106 

140 

Kr-F 

0.26 

0.08 

>  109 

500 

Xenon 

0.175 

0.02 

>  10s 

- 

*Tunable  sources  flamp-driven  ^diode-driven 

YAG  stands  for  Yttrium-Aluminum  Garnet  and  OPO  for  Optical  Parametric 
Oscillator;  rjp  is  pump  laser  efficiency. 
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Formulas 

An  e-m  wave  with  k  ||  B  has  an  index  of  refraction  given  by 
n±  =  [1  -  ujp2e/u)(u  =F  u>ce)]1/2, 

where  ±  refers  to  the  helicity.  The  rate  of  change  of  polarization  angle  6  as  a 
function  of  displacement  s  (Faraday  rotation)  is  given  by 

dO/ds  =  (fc/2)(n_  -  n+)  =  2.36  x  10 4NBf~2  cm-1 , 

where  N  is  the  electron  number  density,  B  is  the  field  strength,  and  /  is  the 
wave  frequency,  all  in  cgs. 

The  quiver  velocity  of  an  electron  in  an  e-m  field  of  angular  frequency  u 
is 

vo  =  eEm&x/muj  =  25.6/1^2Ao  cm  sec-1 

in  terms  of  the  laser  flux  I  =  ci?m2x/87r,  with  I  in  watt/cm2,  laser  wavelength 
Ao  in  yitm.  The  ratio  of  quiver  energy  to  thermal  energy  is 

VFqu/Wth  =  mev02/2kT  =  1.81  x  10_13Ao2//T, 

where  T  is  given  in  eV.  For  example,  if  I  =  1015  W  cm-2 ,  Ao  =  1  ^.m,  T  = 
2keV,  then  VFqu/VFth  ~  0.1. 

Pondermotive  force: 


T  =  N\7(E2)/8irNc, 


where 

Nc  =  1.1  x  1021  A0-2cm-3. 

For  uniform  illumination  of  a  lens  with  /-number  F,  the  diameter  d  at 
focus  (85%  of  the  energy)  and  the  depth  of  focus  l  (distance  to  first  zero  in 
intensity)  are  given  by 

d  fa  2A4F\6/6dl  and  l  fa  ±2F2\9/9dl- 

Here  0  is  the  beam  divergence  containing  85%  of  energy  and  6dl  is  the 
diffraction-limited  divergence: 

9dl  =  2.44A/6, 

where  b  is  the  aperture.  These  formulas  are  modified  for  nonuniform  (such  as 
Gaussian)  illumination  of  the  lens  or  for  pathological  laser  profiles. 
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ATOMIC  PHYSICS  AND  RADIATION 


Energies  and  temperatures  are  in  eV ;  all  other  units  are  cgs  except  where 
noted.  Z  is  the  charge  state  (Z  =  0  refers  to  a  neutral  atom);  the  subscript  e 
labels  electrons.  N  refers  to  number  density,  n  to  principal  quantum  number. 
Asterisk  superscripts  on  level  population  densities  denote  local  thermodynamic 
equilibrium  (LTE)  values.  Thus  Nn*  is  the  LTE  number  density  of  atoms  (or 
ions)  in  level  n. 

Characteristic  atomic  collision  cross  section: 

(1)  7to02  =  8.80  x  10-17cm2. 


Binding  energy  of  outer  electron  in  level  labelled  by  quantum  numbers  n,  l : 

>y2  rp H 

<2>  = 

where  E ^  =  13.6  eV  is  the  hydrogen  ionization  energy  and  A i  =  0.75 1~5 , 
l  5,  is  the  quantum  defect. 

Excitation  and  Decay 

Cross  section  (Bethe  approximation)  for  electron  excitation  by  dipole 
allowed  transition  m  — *•  n  (Refs.  32,  33): 


(3) 


=  2.36  x  10 


—  13  fnm9(.n i  m) 
cA  Enr 


cm 


where  fnrn  is  the  oscillator  strength,  g(n,  m)  is  the  Gaunt  factor,  e  is  the 
incident  electron  energy,  and  A Enm  =  En  —  Em. 

Electron  excitation  rate  averaged  over  Maxwellian  velocity  distribution,  Xmn 
=  Ne(  &mn  V  )  (Refs.  34,  35): 


(4) 


=  1.6  x  10~5 


fnm  (g(n,m))Ne 
A  EnmTl/2 


eXP 


/  AEnm 

V 


-1 
sec  , 


where  (g(n,m))  denotes  the  thermal  averaged  Gaunt  factor  (generally  ~  1  for 
atoms,  ~  0.2  for  ions). 
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Rate  for  electron  collisional  deexcitation: 


(5)  Ynm  =  (Nm*/ Nn*)Xmn. 

Here  iVm* / Nn*  =  (gm / gn)  exp(AEnm /Te)  is  the  Boltzmann  relation  for  level 
population  densities,  where  gn  is  the  statistical  weight  of  level  n. 

Rate  for  spontaneous  decay  n  — >  m  (Einstein  A  coefficient)34 

(6)  Anm  —  4.3  X  10  (9m/5n)/mn(^^'nm)  Sec 


Intensity  emitted  per  unit  volume  from  the  transition  n  ^  ro  in  an  optically 
thin  plasma: 

(7)  Inm  =  1-6  x  10  Ajim Nn AEfim  watt/cm  . 


Condition  for  steady  state  in  a  corona  model: 

(8)  N0Ne(aonv)  =  NnAn0, 

where  the  ground  state  is  labelled  by  a  zero  subscript. 

Hence  for  a  transition  n  —*  m  in  ions,  where  {gin,  0))  «  0.2, 

watt 
cm3 


(9) 


=  5.1  x  10 


-25  /nmflO^e^O 


9m  ^ 


1/2 


/  AEnm  \  3  ( 
V  AEn0  )  6XP  V 


AEn o  \ 
Te  ) 


Ionization  and  Recombination 

In  a  general  time-dependent  situation  the  number  density  of  the  charge 
state  Z  satisfies 

-S(Z)N(Z)  —  a(Z)N(Z) 

+S(Z  -  1)N(Z  -  1)  +  a(Z  +  l)N(Z  +  1)  . 

Here  S(Z)  is  the  ionization  rate.  The  recombination  rate  a(Z)  has  the  form 
oc{Z)  =  ar(Z )  +  Nea3(Z),  where  ar  and  03  are  the  radiative  and  three-body 
recombination  rates,  respectively. 


(10) 


dN(Z ) 
dt 


=  Nf, 
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Classical  ionization  cross-section36  for  any  atomic  shell  j 


(11)  <7i  =  6  x  10  14bjgj(x)/Uj2  cm2. 

Here  bj  is  the  number  of  shell  electrons;  Uj  is  the  binding  energy  of  the  ejected 
electron;  x  =  e/Uj,  where  e  is  the  incident  electron  energy;  and  g  is  a  universal 
function  with  a  minimum  value  <7min  «  0.2  at  x  «  4. 

Ionization  from  ion  ground  state,  averaged  over  Maxwellian  electron  distribu¬ 
tion,  for  0.02  <  Te/E ^  <  100  (Ref.  35): 


(12)  S(Z )  =  10 


-5 


(TJE^)1/2 


(££)3/2(6.0  +  Te/££) 


exp 


(-« 


cm 


'/see, 


where  E ^  is  the  ionization  energy. 

Electron-ion  radiative  recombination  rate  (e  +  N(Z)  — *  N{Z  —  1)  +  hv ) 
for  Te/Z2  <  400  eV  (Ref.  37): 


(13) 


otr(Z)  =  5.2  x  10-14Z 


1/2 


0.43+iln(E^/Te) 


+0.469(E^/Te)"1/3 


3  / 

cm  /sec. 


For  1  eV  <  Te/Z2  <  15 eV,  this  becomes  approximately35 
(14)  ar(Z)  =  2.7  X  10-13Z2Te-1/2cm3/sec. 


OQ 

Collisional  (three-body)  recombination  rate  for  singly  ionized  plasma: 
(15)  c*3  =  8.75  X  10-27Te-4'5  cm6 /sec. 


Photoionization  cross  section  for  ions  in  level  n,  l  (short-wavelength  limit): 

(16)  <Tph(n,  l)  =  1.64  x  10 ~16Z5/n3K7+21  cm2, 

where  K  is  the  wavenumber  in  Rydbergs  (1  Rydberg  =  1.0974  X  105  cm-1). 
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Ionization  Equilibrium  Models 

Saha  equilibrium:39 


(17) 


NeNi^jZ) 

Nn*(Z  -  1) 


=  6.0  x  10 


21  9 1 


Z  rp  3/2 


1 

9n 


exp 


where  g ^  is  the  statistical  weight  for  level  n  of  charge  state  Z  and  £^,(n,  l) 
is  the  ionization  energy  of  the  neutral  atom  initially  in  level  (n,  l),  given  by 
Eq.  (2). 

In  a  steady  state  at  high  electron  density, 


(18) 


NeN*(Z)  S(Z  -  1) 
N*(Z-  1)  “  cT3 


a  function  only  of  T. 

Conditions  for  LTE:39 

(a)  Collisional  and  radiative  excitation  rates  for  a  level  n  must  satisfy 
(19)  Ynrn  ^  lOAnrn. 


(b)  Electron  density  must  satisfy 

(20)  Ne  >  7  x  1018Z7n_17/2(T/Ej>)1/2cm"3. 


Steady  state  condition  in  corona  model: 


(21) 


N(Z  -  1)  _  ar 
N(Z)  ~  S(Z  -  1)' 


Corona  model  is  applicable  if40 

(22)  1012t/_1  <  Ne  <  1016Te7/2  cm-3, 

where  t/  is  the  ionization  time. 
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Radiation 

N.  B.  Energies  and  temperatures  are  in  eV;  all  other  quantities  are  in 
cgs  units  except  where  noted.  Z  is  the  charge  state  (Z  —  0  refers  to  a  neutral 
atom);  the  subscript  e  labels  electrons.  N  is  number  density. 

Average  radiative  decay  rate  of  a  state  with  principal  quantum  number  n  is 
(23)  An  =  ^2  Anm  =  1.6  X  10 10Z4n_9/2  sec. 

m<n 


Natural  linewidth  (A E  in  eV): 

(24)  AE  At  =  h  =  4.14  x  10“15eVsec, 

where  At  is  the  lifetime  of  the  line. 

Doppler  width: 

(25)  AA/A  =  7.7  x  10 _5(T/p)1/2, 

where  p  is  the  mass  of  the  emitting  atom  or  ion  scaled  by  the  proton  mass. 
Optical  depth  for  a  Doppler-broadened  line:39 

(26)  r  =  3.52 x  10~13 fnrn\(Mc2 /kT)1^2 NL  =  5.4x10 “9/mn A(p/T)1/2iVL, 

where  fnm  is  the  absorption  oscillator  strength,  A  is  the  wavelength,  and  L  is 
the  physical  depth  of  the  plasma;  M,  N ,  and  T  are  the  mass,  number  density, 
and  temperature  of  the  absorber;  p  is  M  divided  by  the  proton  mass.  Optically 
thin  means  r  <  1. 

Resonance  absorption  cross  section  at  center  of  line: 

(27)  <r\=\c  =  5.6  x  10_13A2/AA cm2. 

Wien  displacement  law  (wavelength  of  maximum  black-body  emission): 

(28)  Amax  =  2.50  X  10-5T-1  cm. 

Radiation  from  the  surface  of  a  black  body  at  temperature  T: 

(29)  W  =  1.03  x  105T4  watt/cm2. 
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Bremsstrahlung  from  hydrogen-like  plasma:26 

(30)  PBr  =  1-69  x  10 ~32NeTe1/2  [ Z2N(Z )]  watt/cm3, 


where  the  sum  is  over  all  ionization  states  Z. 
Bremsstrahlung  optical  depth:41 


(31)  t  =  5.0  X  10~38  NeNiZ2gLT~7/2, 

where  g  «  1.2  is  an  average  Gaunt  factor  and  L  is  the  physical  path  length. 

Inverse  bremsstrahlung  absorption  coefficient42  for  radiation  of  angular  fre¬ 
quency  uj: 


(32)  k  =  3.1  X  10  7Zne2\nAT  372cu  2(1  —  uj2/u>2)  172  cm  1; 

here  A  is  the  electron  thermal  velocity  divided  by  V,  where  V  is  the  larger  of 
uj  and  Up  multiplied  by  the  larger  of  Ze2/kT  and  h/  (mkT)1^2 . 

Recombination  (free-bound)  radiation: 


(33) 


Pr  =  1.69  X  10~ 32 NeT( 


1/2 


£  [*•»<»  (+) 


watt/cm" 


Cyclotron  radiation26  in  magnetic  field  B: 

(34)  Pc  =  6.21  x  10 ~28B2NeTe  watt/cm3. 

For  NekTe  =  NikTi  =  B2 /16tc  (/3  =  1,  isothermal  plasma),26 

(35)  Pc  =  5.00  X  10 ~38N2T2  watt/cm3. 


Cyclotron  radiation  energy  loss  e-folding  time  for  a  single  electron:41 

9.0  x  108P-2 


(36) 


2.5  +  7 


sec, 


where  7  is  the  kinetic  plus  rest  energy  divided  by  the  rest  energy  me2. 
Number  of  cyclotron  harmonics41  trapped  in  a  medium  of  finite  depth  L: 

(37)  mtr  =  (57  (3BL)1/6, 

where  (3  =  SirNkT/B2. 

Line  radiation  is  given  by  summing  Eq.  (9)  over  all  species  in  the  plasma. 
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ATOMIC  SPECTROSCOPY 


Spectroscopic  notation  combines  observational  and  theoretical  elements. 
Observationally,  spectral  lines  are  grouped  in  series  with  line  spacings  which 
decrease  toward  the  series  limit.  Every  line  can  be  related  theoretically  to  a 
transition  between  two  atomic  states,  each  identified  by  its  quantum  numbers. 

Ionization  levels  are  indicated  by  roman  numerals.  Thus  C I  is  unionized 
carbon,  CII  is  singly  ionized,  etc.  The  state  of  a  one-electron  atom  (hydrogen) 
or  ion  (Hell,  Li  III,  etc.)  is  specified  by  identifying  the  principal  quantum 
number  n  =  1,2,...,  the  orbital  angular  momentum  1  =  0, 1, . . . ,  n  —  1,  and 
the  spin  angular  momentum  s  =  ±^.  The  total  angular  momentum  j  is  the 
magnitude  of  the  vector  sum  of  1  and  s,  j  =  l  ±  \  ( j  >  ^).  The  letters  s, 
p,  d,  f,  g,  h,  i,  k,  1,  .  . . ,  respectively,  are  associated  with  angular  momenta 
1  =  0,  1,  2,  3,  4,  5,  6,  7,  8,  ....  The  atomic  states  of  hydrogen  and  hydrogenic 
ions  are  degenerate:  neglecting  fine  structure,  their  energies  depend  only  on  n 
according  to 

R00hcZ2n~ 2  RyZ2 

n  1  +  m/M  n2  ’ 

where  h  is  Planck’s  constant,  c  is  the  velocity  of  light,  m  is  the  electron  mass, 
M  and  Z  are  the  mass  and  charge  state  of  the  nucleus,  and 

Roo  =  109,  737cm-1 

is  the  Rydberg  constant.  If  En  is  divided  by  he,  the  result  is  in  wavenumber 
units.  The  energy  associated  with  a  transition  m  —*  n  is  given  by 

A Emn  -  Ry(l/m2  -  1/n2), 

with  m  <  n  (m  >  n)  for  absorption  (emission)  lines. 

For  hydrogen  and  hydrogenic  ions  the  series  of  lines  belonging  to  the 
transitions  m  — >  n  have  conventional  names: 


Transition 

1  — *  n 

2  — ►  n 

C 

T 

CO 

4  — ►  n 

5  — +  n 

6  —*■  n 

Name 

Lyman 

Balmer 

Paschen 

Brackett 

Pfund 

Humphreys 

Successive  lines  in  any  series  are  denoted  a,  (3,  7,  etc.  Thus  the  transition  1  — > 
3  gives  rise  to  the  Lyman-/?  line.  Relativistic  effects,  quantum  electro  dynamic 
effects  (e.g.,  the  Lamb  shift),  and  interactions  between  the  nuclear  magnetic 
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moment  and  the  magnetic  field  due  to  the  electron  produce  small  shifts  and 
splittings,  ^  10-2cm-1;  these  last  are  called  “hyperfine  structure.” 

In  many-electron  atoms  the  electrons  are  grouped  in  closed  and  open 
shells,  with  spectroscopic  properties  determined  mainly  by  the  outer  shell. 
Shell  energies  depend  primarily  on  n;  the  shells  corresponding  to  n  =  1,  2, 
3, . . .  are  called  K ,  L,  M,  etc.  A  shell  is  made  up  of  subshells  of  different 
angular  momenta,  each  labeled  according  to  the  values  of  n,  l,  and  the  number 
of  electrons  it  contains  out  of  the  maximum  possible  number,  2(2 1  +  1).  For 
example,  2p5  indicates  that  there  are  5  electrons  in  the  subshell  corresponding 
to  l  =  1  (denoted  by  p)  and  n  —  2. 

In  the  lighter  elements  the  electrons  fill  up  subshells  within  each  shell 
in  the  order  s,  p,  d,  etc.,  and  no  shell  acquires  electrons  until  the  lower  shells 
are  full.  In  the  heavier  elements  this  rule  does  not  always  hold.  But  if  a 
particular  subshell  is  filled  in  a  noble  gas,  then  the  same  subshell  is  filled  in 
the  atoms  of  all  elements  that  come  later  in  the  periodic  table.  The  ground 

state  configurations  of  the  noble  gases  are  as  follows: 

He  Is2 

Ne  ls22s22p6 

Ar  ls22s22p63s23p6 

Kr  Is2  2s2  2p6  3s2  3p6  3d104s2  4p6 

Xe  ls22s22p63s23p63d104s24p64d105s25p6 

Rn  ls22s22p63s23p63d104s24p64d104f145s25p65d106s26p6 


Alkali  metals  (Li,  Na,  K,  etc.)  resemble  hydrogen;  their  transitions  are  de¬ 
scribed  by  giving  n  and  l  in  the  initial  and  final  states  for  the  single  outer 
(valence)  electron. 

For  general  transitions  in  most  atoms  the  atomic  states  are  specified  in 
terms  of  the  parity  (  — l)Eii  and  the  magnitudes  of  the  orbital  angular  momen¬ 
tum  L  =  Xlj,  the  spin  S  =  Xs*,  and  the  total  angular  momentum  J  =  L  +  S, 
where  all  sums  are  carried  out  over  the  unfilled  subshells  (the  filled  ones  sum 
to  zero).  If  a  magnetic  field  is  present  the  projections  Ml,  Ms,  and  M  of 
L,  S,  and  J  along  the  field  are  also  needed.  The  quantum  numbers  satisfy 
\Ml\  <  L  <  vl,  |Ms|  <  S  <  v/2,  and  \M\  <  J  <  L  +  S,  where  v  is  the 
number  of  electrons  in  the  unfilled  subshell.  Upper-case  letters  S,  P,  D,  etc., 
stand  for  L  =  0,1,  2,  etc.,  in  analogy  with  the  notation  for  a  single  electron. 
For  example,  the  ground  state  of  Cl  is  described  by  3p5  2Pg/2-  The  first  part 
indicates  that  there  are  5  electrons  in  the  subshell  corresponding  to  n  =  3  and 
1  =  1.  (The  closed  inner  subshells  Is2 2s2  2p6 3s2,  identical  with  the  configura¬ 
tion  of  Mg,  are  usually  omitted.)  The  symbol  ‘P’  indicates  that  the  angular 
momenta  of  the  outer  electrons  combine  to  give  L  =  1.  The  prefix  ‘2’  repre¬ 
sents  the  value  of  the  multiplicity  2 S  +  1  (the  number  of  states  with  nearly  the 
same  energy),  which  is  equivalent  to  specifying  S  =  5.  The  subscript  3/2  is 
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the  value  of  J.  The  superscript  ‘o’  indicates  that  the  state  has  odd  parity;  it 
would  be  omitted  if  the  state  were  even. 

The  notation  for  excited  states  is  similar.  For  example,  helium  has  a  state 
ls2s  3Si  which  lies  19.72  eV  (159,856  cm-1)  above  the  ground  state  Is2  1So- 
But  the  two  “terms”  do  not  “combine”  (transitions  between  them  do  not  occur) 
because  this  would  violate,  e.g.,  the  quantum-mechanical  selection  rule  that 
the  parity  must  change  from  odd  to  even  or  from  even  to  odd.  For  electric 
dipole  transitions  (the  only  ones  possible  in  the  long- wavelength  limit),  other 
selection  rules  are  that  the  value  of  l  of  only  one  electron  can  change,  and  only 
by  A  l  =  ±1;  AS  =  0;  A  L  =  ±1  or  0;  and  A  J  =  ±1  or  0  (but  L  =  0  does  not 
combine  with  L  =  0  and  J  =  0  does  not  combine  with  J  —  0).  Transitions 
are  possible  between  the  helium  ground  state  (which  has  S  =  0,  L  =  0,  J  =  0, 
and  even  parity)  and,  e.g.,  the  state  ls2p  1P°  (with  S  =  0,  L  =  1,  J  =  1, 
odd  parity,  excitation  energy  21.22 eV).  These  rules  hold  accurately  only  for 
light  atoms  in  the  absence  of  strong  electric  or  magnetic  fields.  Transitions 
that  obey  the  selection  rules  are  called  “allowed” ;  those  that  do  not  are  called 
“forbidden.” 

The  amount  of  information  needed  to  adequately  characterize  a  state  in¬ 
creases  with  the  number  of  electrons;  this  is  reflected  in  the  notation.  Thus43 
O II  has  an  allowed  transition  between  the  states  2p  3p 
2f7°/2  and  2p2(1D)3d/  2F7/2  (and  between  the  states  obtained  by  changing 

J  from  7/2  to  5/2  in  either  or  both  terms).  Here  both  states  have  two  elec¬ 
trons  with  n  =  2  and  1  =  1;  the  closed  subshells  Is2 2s2  are  not  shown.  The 
outer  (n  =  3)  electron  has  l  =  1  in  the  first  state  and  l  =  2  in  the  second. 
The  prime  indicates  that  if  the  outermost  electron  were  removed  by  ionization, 
the  resulting  ion  would  not  be  in  its  lowest  energy  state.  The  expression  (XD) 
give  the  multiplicity  and  total  angular  momentum  of  the  “parent”  term,  i.e., 
the  subshell  immediately  below  the  valence  subshell;  this  is  understood  to  be 
the  same  in  both  states.  (Grandparents,  etc.,  sometimes  have  to  be  specified 
in  heavier  atoms  and  ions.)  Another  example43  is  the  allowed  transition  from 
2p2(3P)3p  2Pi/2  (or  2p3/2)  to  2p2(1D)3d'  2Si/2,  in  which  there  is  a  “spin 

flip”  (from  antiparallel  to  parallel)  in  the  n  =  2,  l  =  1  subshell,  as  well  as 
changes  from  one  state  to  the  other  in  the  value  of  l  for  the  valence  electron 
and  in  L. 

The  description  of  fine  structure,  Stark  and  Zeeman  effects,  spectra  of 
highly  ionized  or  heavy  atoms,  etc.,  is  more  complicated.  The  most  important 
difference  between  optical  and  X-ray  spectra  is  that  the  latter  involve  energy 
changes  of  the  inner  electrons  rather  than  the  outer  ones;  often  several  electrons 
participate. 
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AFTERWORD 


The  NRL  Plasma  Formulary  originated  over  twenty  years  ago  and 
has  been  revised  several  times  during  this  period.  The  guiding  spirit  and  person 
primarily  responsible  for  its  existence  and  upkeep  is  Dr.  David  Book.  The 
Formulary  has  been  set  in  T^K  by  Dave  Book,  Todd  Brun,  and  Robert  Scott. 
I  am  indebted  to  Dave  for  providing  me  with  the  TgX  files  for  the  Formulary 
and  his  assistance  in  its  re-issuance.  Finally,  I  thank  readers  for  communicating 
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